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I.  INTRODUCTION 


During  the  period  of  this  program  the  scientific  objectives  were  to  complete  a  study 
of  one-dimensional  electrostatic  plasmas  using  the  discrete  Hamiltonian  method,  amd 
to  construct  a  model  of  ridged  high  harmonic  gyrotron  oscillators  using  '‘standard^ 
gyrotron  modeling  techniques.  Ridged  gyrotrons  appear  ideal  for  the  future  applica¬ 
tion  of  the  multiple  time  scale,  discrete  Hamiltonian  method.  As  with  all  gyrotrons, 
ridged  gyrotrons  are  characterized  by  several  primary  time  scales,  which  are  the  cy¬ 
clotron  period,  the  cavity  RF  field  oscillation  period,  and  the  transit  time  through  the 
gyrotron  tube.  If  only  one  harmonic  is  being  considered,  then  the  cyclotron  period  and 
the  cavity  oscillation  period  need  not  be  treated  separately  since  only  the  difference 
between  them  has  significance.  For  high  harmonic  emission,  however,  the  ridged  gy¬ 
rotron  appears  to  require  the  simultaneous  treatment  of  multiple  cyclotron  harmonics 
and  hence  multiple  harmonic  time  scales.  Jn  addition,  mode  competition,  which  can 
lead  to  coupling  between  multiple  RF,.cavity  modes,  is  possible  and  would  also  result 

t  \ 

in  the  need  to  treat  several  cavity  frequencies  simultaneously.  *'  '  J  • 

In  investigating  electrostatic  plasmas,  difference  equations  were  obtained  by  inte¬ 
grating  over  unperturbed  orbits,  allow  time  steps  much  longer  than  the  linear  mode 
oscillation  period.  In  our  gyrotron  modeling,  we  developed  a  detailed  lineu  theory,  and 
a  nonlinear  model  using  a  non-interacting  test  pairticle  approach  to  solving  the  electron 
guiding  center  coordinate  system  equations  of  motion  in  assumed  cold  cavity  RF  fields 
for  ridged  cavities.  The  ridged  cavity  systems  which  we  studied  have  good  potential  for 
producing  high  frequency  radiation  from  high  harmonic  emission.  Our  studies  can  be 
used  to  provide  direct  theoretical  support  for  the  ridged  rectangulu  cavity  experimen¬ 
tal  program  being  conducted  by  Ferendeci  [l].  Two  papers  were  submitted  for  journal 
publication  during  this  contract  period  and  one  paper  was  presented  at  the  Eleventh 
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International  Conference  on  Infrared  and  Millimeter  Waves.  Copies  of  the  papers  are 
given  in  the  Appendices. 


n.  TECHNICAL  REPORT 


A.  DISCRETE  HAMILTONIAN  METHOD 

In  order  to  test  the  discrete  Hamiltonian  method  in  relatively  simple  circumstances^ 
Dr.  Decyk  and  Dr.  Menyuk  studied  one-dimensional  electrostatic  plasmas,  which  were 
either  Maxwellian  distributed  or  had  a  warm  beam  present.  The  discrete  Hamiltonian 
method  is  predicated  on  two  basic  ideas:  1)  By  integrating  over  unperturbed  orbits,  it  is 
possible  to  take  time  steps  which  are  substantially  longer  than  is  possible  using  leapfrog 
codes  which  simply  extrapolate  linearly  in  phase  space.  2)  The  equations  of  motion 
of  a  closed  plasma  (including  electromagnetic  interactions)  are  Hamiltonian.  There¬ 
fore,  by  retaining  this  property  in  the  discrete  equations,  one  avoids  certain  nonlinear 
instabilities  which  can  be  very  damaging  to  the  long-term  behavior  of  !^he  system. 

In  principle,  the  discrete  Hamiltonian  equations  can  be  written  such  that  the 
discrete  change  in  the  particle  spatial  coordinate  and  momentum  Az  and  Ap  depend 
on  the  initial  value  of  x  amd  the  final  value  of  p.  When  integrating  over  the  unperturbed 
orbits,  we  do  not  assume  that  the  field  is  frozen,  the  assumption  which  is  generally 
made  in  leapfrog  codes.  Instead,  we  assume  that  the  waves  vary  according  to  the  linear 
dispersion  relation.  This  assumption  allows  us  in  principle  to  take  time  steps  which 
are  long  compared  to  the  linear  mode  oscillation  period,  but  are  short  compared  to 
the  nonlinear  time  scale.  Our  formulation  requires  us  to  sum  over  the  different  linear 
modes  in  the  system.  In  the  one-dimensional,  electrostatic  problem  which  we  have 
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studied,  the  discrete  equations  describing  the  particle  push  may  be  written  as  [2] 


_  t;**  =  —  ^ 


^  m(v"+i  -wjk/*) 
-  cos  (hz"  + 


I  cos  [ib**  +  -  £4;*)  r  + 


-ft+l 


-  »"  =  „»+>r  -  •£  .in  r  +  «] 

—  sin  {kz^  +  <f>l)  —  —  w*)  r  cos  [fcz"  +  —  w*)  r  +  ^jj]  | , 


(1) 


where  the  superscript  n  indicates  the  time  step,  the  subscript  k  indicates  the  mode 
number,  and  the  ordinary  variables  k  and  u  are  the  wavenumber  and  frequency  respec¬ 
tively.  The  variables  r,  ^!ci  respectively  the  size  of  the  time  step,  the 

mode  potential,  the  mode  phase,  and  the  phase  velocity.  In  the  sum  over  k,  we  must 
distinguish  between  forward-going  and  backward-going  waves.  We  do  so  by  use  of  the 
particle  current  as  well  as  the  charge.  Writing 


=  (A;*/45r)  cos  (-Wifcr  +  4>±k) , 
=  {k^/4ir)  $±fcsm(-w±jkr  +  <f>±k) , 


(2) 


we  find 

—  {^/^)[Pkr ^jkrl^k)7  ^  =  (l/2)(Pfcr  ~  ^J*r/Wfc), 

(3) 

=  (l/2](pj|,j  —  kjici/(*ik)7  S  =  {l/2){pki  +  fcjfci/wfc), 

where  -h  and  —  indicate  respectively  the  forward-going  and  backward-going  waves  while 
the  subscripts  r  and  i  indicate  the  real  and  imaginary  parts  of  the  chairge  auid  current 
in  the  usual  Fourier  decomposition. 

Equations  (1-3)  are  only  weakly  implicit;  that  is  to  say  that  the  discrete  equation 
depends  on  the  field  at  the  beginning  of  the  time  step.  As  a  consequence,  one  can  solve 
for  the  change  in  z  and  v  directly  for  each  particle  using  Newton’s  method. 


I 
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It  is  natural  to  inquire  whether  the  lack  of  time-centering  of  the  fields  which  this 
represents  has  any  bad  effects.  It  is  possible,  in  fact  to  vrrite  down  a  strongly  implicit 
form  of  (1)  which  is  time-centered  and  strongly  implicit 


jfi+l  _gn  _  yn+lj. 


-  -  Ult)  Tcoa  [ks*  +  -  w*)  r  + 


(4) 

The  linear  dispersion  relation  which  results  from  (4)  agrees  exactly  with  the  linear 
dispersion  relation  obtained  by  solving  the  continuous,  exact  equations  of  motion,  in 
contrast  to  (1)  which  has  small  differences.  We  tried  solving  (4)  using  a  predictor- 
corrector  approach.  Remarkably,  however,  this  approach  did  not  converge.  We  will 
return  to  the  reason  for  this  later.  As  a  consequence  of  this  non-convergence,  we 
always  used  (1)  in  practice. 

The  first  problem  we  faced  in  solving  these  equations  is  that  calculating  all  of 
the  sine  and  cosine  functions  in  (1)  for  each  of  the  modes  proved  to  be  prohibitively 
expensive  even  for  a  restricted  number  of  modes  (16  or  32  in  practice).  To  deal  with 
this  difficulty,  we  wrote  the  code,  including  the  Newton’s  method  iteration,  to  take 
advantage  of  vectorization  on  the  CRAY  computers  (where  we  ran  our  code)  and  also 
to  make  maximum  use  of  recursion  relations.  Through  the  use  of  these  relations,  it  is 
possible  to  reduce  the  number  of  cosine  and  sine  evaluations  to  one  each  per  particle  per 
time  step.  With  these  changes,  each  time  step  in  our  code  was  a  factor  of  5  longer  than 
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in  an  optimized  leapfrog  code.  Thus,  the  crossover  point  to  obtain  computational  time 
saving  occurs  when  si  1,  a  value  which  we  regularly  exceeded  in  our  calculations. 

Menyuk  and  Becyk  [3]  have  already  reported  that  when  modes  are  present  with 
phase  velocities  in  the  bulk  of  the  plasma,  the  discrete  Hamiltonian  approach  is  strongly 
unstable.  The  reason  is  that  the  wave  amplitudes  fluctuate  rapidly  on  a  time  scale 
about  equal  to  Up,  and  the  code  cannot  follow  these  fluctuations  since  they  violate  the 
assumption  that  wave  amplitudes  change  slowly  on  the  linear  time  scale.  While  the 
fluctuations  are  due  to  the  use  of  discrete  particles,  little  improvement  is  obtained  in  a 
thermal  simulation  plasma  by  increasing  the  number  of  particles.  The  time  scale  of  the 
fluctuations  scales  like  the  trapping  times  in  the  individual  k-modes,  T/iuet  oc  ■> 

and  Ek  «  where  N  is  the  number  of  particles.  Hence,  r//uci  oc  and 

increases  only  weakly  with  N.  Our  constraint  is  r  <  T/iuct  for  stability.  We  have 
verified  that  increasing  the  number  of  particles  in  our  simulation  by  a  factor  16  allows 
approximately  double  the  size  of  r  before  the  strong  instability  sets  in.  When  it  is 
present,  it  occurs  on  a  time  scale  (jJptT  C2l0-20  for  the  cases  we  ran,  where  T  is  the 
total  rim  time. 

During  this  contract  period,  we  also  found  a  new  source  of  numerical  instability 
which  occurs  on  a  very  long  time  scale,  u/peT  300  when  UpgT  =  4  and  UpeT  400 
when  UpeT  =  3.  This  instability  and  its  effect  on  phase  space  is  shown  in  Figure  1 
in  the  case  Wpgf  =  3.  We  at  first  thought  that  this  problem  was  due  to  the  lack  of 
time  centering  in  the  fields,  but,  as  previously  noted,  taking  predictor-corrector  steps 
to  center  the  fields  did  not  converge  and  thus  only  exacerbated  the  problem.  The  real 
culprit  turned  out  to  be  rapid  fluctuations  in  the  field  amplitudes  due  to  the  random 
particle  assortment  which  still  occurs  even  when  the  waves’  phase  velocities  are  outside 
the  bulk  plasma.  Indeed,  these  rapid  fluctuations  also  account  for  the  failure  of  the 
predictor-corrector  scheme  to  converge.  To  deal  with  this  difiiculty  effectively,  it  is 
necessary  to  use  quiet  start  techniques  [4].  The  dramatic  improvement  which  results 
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can  be  seen  by  comparing  Figure  2  with  Figure  1.  It  should  be  noted,  however,  that 
the  use  of  quiet  start  techniques  is  not  a  solution  over  indefinitely  long  times.  At 
ijJpeT  =  450,  the  particles  have  clearly  been  randomized.  On  a  longer  time  scale,  we 
anticipate  that  the  plasma  will  once  more  be  numerically  unstable.  This  problem  can 
in  principle  be  eliminated  by  periodic  regridding  [4].  We  have  not,  however,  explored 
this  question  ourselves,  feeling  that  this  instability  occurs  on  such  a  long  time  scale 
that  it  is  unlikely  to  be  of  interest  in  most  czises. 

Finally,  we  have  studied  the  beam-plasma  instability,  using  a  10%,  a  1%,  and  a 
0.1%  beam.  The  beam  amplitude  is  lowered  by  decreasing  q  and  m  for  each  particle 
in  the  beam,  holding  q/m  constant,  rather  that  by  decreasing  the  number  of  puticles, 
in  order  to  obtain  a  reasonable  beam  resolution.  The  phase  space  evolution  is  shown 
in  Figure  3  for  the  10%  beam.  The  trapping  and  thermalization  of  the  beam  is  readily 
apparent.  In  Figure  4,  we  show  the  evolution  of  the  field  energy  for  the  three  different 
amplitude  beams.  As  expected,  the  instability  becomes  less  energetic  when  going  from  a 
10%  beam  to  a  1%  beam,  and  disappears  for  a  0.1%  beam.  These  results  are  consistent 
with  the  results  of  a  comparison  leapfrog  code  which  we  have  run. 

We  conclude  that  while  the  discrete  Hamiltonian  method  is  not  simple  to  imple¬ 
ment  in  practice,  it  can  be  useful  in  situations  where  the  number  of  modes  needed  is 
small  2uid  long  time  steps  are  required. 
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B.  DEVELOPMENT  OF  GYROTRON  MODELING  TOOLS 


The  first  stage  in  the  implementation  of  the  discrete  Hamiltonian  method  to  gy- 
rotron  modeling  is  the  choice  of  a  device  to  study  and  the  development  of  a  set  of  basic 
gyrotron  models,  based  upon  conventional  approocimations  and  techniques,  which  can 
be  modified  and  expanded  upon.  As  initial  devices  to  investigate,  we  have  chosen  a 

1  high  harmonic  ridged  rectangular  cavity  and  a  ridged  co-axial  cylindrical  cavity  gy- 

f 

rotron  oscillator.  The  cavity  cross  sections  for  the  two  devices  are  shown  in  Figures 
5a  and  5b.  The  emission  harmonics  are  defined  by  sQ/'i  =  u,  where  u  is  the  cavity 
frequency  and  s  is  the  harmonic  number.  As  emission  for  the  gyrotron  occurs  at  har¬ 
monics  of  the  electron  cyclotron  frequency,  high  frequency  emission  implies  either  the 
use  of  large  magnetic  fields  or  high  harmonic  operation.  If  frequencies  as  high  as  the 
millimeter  wave  regime  are  to  be  achieved  using  conventional  magnets,  one  must  use 
third  and  higher  harmonics.  Studies  of  gyrotrons  with  ridged  rectangular  tubes  [1,  5-9] 
and  ridged  cylindrical  tubes  [10-18]  indicate  that  the  ridged  cavity  structure  provides 
an  excellent  high  harmonic  interaction  and  good  mode  selection. 

To  aid  us  in  our  understanding  of  ridged  gyrotron  oscillators,  we  have  developed 
a  detadled  linear  theory  for  the  interaction  between  the  electron  beam  and  cavity  RF 
field.  Our  model  treats  gyrotrons  with  ridged  cross-sections  whose  basic  symmetry  is 
either  rectangular  or  cylindrical.  The  model  calculates  the  start-oscillation  condition, 
QPb,  where  Q  is  the  cavity  <)uality  factor  and  is  the  initial  beam  power.  The 
start-oscillation  condition  determines  the  beam  power  at  which  a  cavity  will  oscillate 
and  is  needed  both  for  an  oscillator  design  and  to  investigate  mode  competition.  The 
linear  model  idso  solves  for  the  shift  in  the  cavity  resonance  frequency  u)  caused  by 
the  presence  of  the  electron  beam.  It  provides  an  excellent  basis  for  benchmarking  our 
nonlinear  codes,  and  is  also  an  aid  in  the  design  of  an  optimum  high  harmonic  ridged 
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cavity  gyrotron.  Details  of  the  model  are  given  in  the  Appendices  which  contains  papers 
i  which  have  been  submitted  for  journal  publication. 

With  our  linear  model  in  hand,  we  turned  to  the  ridged  rectangular  gyro-TWT 
amplifier  work  of  Ferendeci  and  Han  [1,  5-8]  and  the  ridged  cylindrical  gyrotron  os- 
^  dilator  modeling  of  Chu  and  Dialetis  [14-15]  to  begin  a  program  of  investigation  of 

ridged  cavities.  A  comparison  of  our  modeling  of  the  ridged  cavity  RF  field  with  that 
of  Ferendeci  and  Han  revealed  that  while  there  was  exact  agreement  between  our  ap- 
I  proach  and  the  basic  equations  describing  the  RF  field,  the  solutions  for  several  of  the 

resonant  frequencies  differed.  Upon  a  careful  examination  of  the  RF  field  described  in 
Han’s  PhD  thesis  [5],  we  have  concluded  that  extraneous  terms  had  been  included  by 
Han  in  the  solution  of  the  cavity  frequency  dispersion  relation.  For  the  TE02  mode 
which  Ferendeci  and  Han  believed  to  be  optimal  for  their  design,  no  discrepancies  exist 
with  our  calculations. 

We  find,  in  agreement  with  Ferendeci  and  Han,  that  the  addition  of  ridges  greatly 
enl.ances  the  interaction  between  the  electron  beam  and  the  RF  field  at  high  harmonics. 
Figure  6  shows  the  TE021  start  oscillation  condition  as  a  function  of  magnetic  field  Bo 
for  a  smooth  walled  rectangular  cavity  suggested  by  Ferendeci  and  Han.  See  Figure 
5b.  The  electron  beam  is  placed  in  the  center  of  the  cavity  to  allow  as  large  a  Larmor 
radius  as  possible.  For  the  smooth  walled  cavity  and  the  chosen  beam  position,  there 
is  emission  only  at  even  harmonics.  With  increasing  harmonic  number  the  gyrotron 
interaction  rapidly  weakens,  and  the  required  start  oscillation  condition  increases.  The 
addition  of  ridges  to  the  cavity  creates  strong  fringe  fields  about  the  ridge  openings  as 
is  evident  from  Figure  7.  These  fringe  fields  lead  to  a  strong  high  harmonic  emission.  In 
Figure  8,  we  show  the  result  that  adding  ridges  has  upon  the  start  oscillation  condition. 
The  gyrotron  emission  interaction  now  occurs  at  all  cyclotron  harmonics.  The  high 
harmonic  start  oscillation  condition  is  very  low,  and  even  decreases  with  increasing 
harmonic  numbw.  Ferendeci  and  Han  [1,  5-8]  did  not  analyse  their  equations  and 
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did  not  find  this  effect.  Hannonics  greater  than  6  for  the  TE^oai  mode  have  Laxmor 
radii  which  would  cause  the  beam  electrons  to  strike  the  cavity  wall  and  thus  are  not 
included  in  the  figure. 
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Modes  other  than  the  TEoai  mode  may  be  excited  by  the  electron  beam.  In  Figure 
9  we  show  the  start  oscillation  condition  for  the  modes  TEmni  where  m  =  0, 1, 2  and 
n  =  1, 2, 3, 4, 5.  What  is  evident  from  this  figure  is  that  other  modes  have  high  harmonic 
resonances  with  start  oecillation  conditions  smaller  that  that  of  the  sixth  harmonic  of 
TEoai  •  Also,  the  harmonic  emission  resonances  of  TEoai  ere  likely  to  suffer  from  mode 
competition  from  the  other  strong  emission  resonances.  As  our  treatment  of  the  RF 
cavity  field  modes  is  approximate  (only  one  expansion  term  is  used  in  the  ridges)  the 
cavity  frequencies  for  the  modes  and  therefore  the  values  of  the  background  magnetic 
field  for  the  harmonic  resonances  cannot  be  exactly  determined  from  our  present  theory. 
The  relative  positions  of  the  harmonics  in  Figure  9  shovid  therefore  be  taken  as  an 
indication  of  possible  mode  competition  and  not  a  definitive  statement  as  to  where  it 
will  exist. 

Start  oscillation  calculations  were  conducted  for  the  ridged  cylindrical  cavity  as 
well.  In  Figure  10  we  show  for  an  axis  encircling  beam  how  the  presence  of  the  ridges 
can  greatly  reduce  the  value  of  QP^.  We  found  that  our  model  agrees  with  the  value 
of  QPi  calculated  for  the  electron  cyclotron  maser  interaction  by  Chu  and  Dialetis 
[14-15].  However,  Chu  and  Dialetis  did  not  include  the  peniotron  resonance  in  their 
modeling.  The  peniotron  interaction,  which  can  lead  to  a  very  high  peak  efficiency 
[19],  is  always  present  in  gyrotron  devices  along  with  the  standard  electron  cyclotron 
maser  interaction.  As  the  peniotron  interaction  was  actually  the  stronger  interaction 
over  much  of  the  parameter  space  considered  by  Chu  and  Dialetis,  the  values  given  in 
their  paper  for  the  minimum  start  oscillation  condition  are  not  correct. 

In  addition  to  the  linear,  ridged  model,  we  have  also  developed  a  nonlinear,  ridged 
gyrotron  model.  The  nonlinear  model  uses  a  non-interacting  test  particle  approach  to 
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solving  the  electron  equation  of  motion 


dz 


(5) 


using  the  assumed  cold  cavity  RF  fields.  (See  the  Appendices  for  a  description  of 
the  ridged  cavity  fields.)  In  (5)  is  the  dimensionless  momentum  vector,  U  =  7/9, 
where  7  give  the  particle  energy  and  0  is  the  velocity  in  terms  of  the  speed  of  light. 
Dimensionless  tmits  scaled  to  the  cavity  length  R  are  used  in  the  modeling  so  that 
the  results  are  able  to  be  scaled  to  any  frequency  cavity.  The  nonlinear  code  uses  the 
standard  gyrotron  modeling  approach  in  that  the  RF  field  in  expanded  in  terms  of  the 
cyclotron  harmonics  in  an  infinite  sine-Bessel  series  about  the  electron  guiding  center. 
In  the  code,  this  sum  is  truncated  to  a  small  set  of  terms,  typically  3-10.  If  only  one 
harmonic  expansion  term  is  dominant,  then  using  the  fact  that  (a;  —  8Cl/'))t  C  ujt,  a 
slow  time  scale  treatment  can  be  used  to  take  time  steps  much  larger  that  the  cyclotron 
period.  Comparison  in  the  weak  field  limit  has  shown  good  agreement  between  the 
linear  theory  and  the  nonlinear  model. 
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rV.  FIGURES 


Figure  1.  Development  of  the  weak  numerical  instability  due  to  wave  amplitude  fluctuations 
{upgT  =  Z,  N  =  40,000,  rimode  —  16)*  *)  =  300,  b)  Up^T  =  450,  c)  field 

energy  vs.  time. 

Figure  2.  Effect  of  a  quiet  start  on  the  weak  numerical  instability  (same  parameters  as  in 
Figure  1).  a)  u>ptT  =  0,  b)  =  450,  c)  field  energy  vs.  time.  Note  at  T  =  0, 
the  regular  pattern  characteristic  of  quiet  start  simulations. 

Figure  3.  Phase  space  plots  of  the  beam-plasma  instability  (10%  beam,  Ntotai  =  32,000, 
Nheam  —  3,200,  <itpe^  —  2,  flmode  —  1®)*  *)  ~  0*  b)  itfptT  =  50,  c)  Up^T  =  100. 

Figure  4.  Field  energy  vs.  time,  a]  10%  beam,  b)  1%  beam,  c)  0.1%  beam. 

Figure  5a.  Cross-section  of  the  ridged  cylindrical  gyrotron. 

Figure  5b.  Cross-section  of  the  ridged  rectangular  gyrotron. 

Figure  6.  TEoai  mode  start  oscillation  condition  for  a  smooth  walled  rectangular  cavity  for 
an  electron  beam  with  energy  72  keV  and  Pg  =  0.28.  The  cavity  parameters  are 
Lx  =  2.1987i2  and  L,  =  15i2.  The  resonwces  are  labeled  by  their  hMmonic 
number. 

Figure  7.  Eg  RF  field  for  a  ridged  rectangular  cavity  with  L*  =  2.1987i?,  Ro  =  1.1983/2, 
V)  =  0.05455/2,  iadN  =  3. 

Figure  8.  TEoji  mode  start  oscillation  condition  for  a  the  ridged  rectangulair  cavity  for  an 
electron  beam  with  72  keV  and  0,  =  0.28.  The  resonances  are  labeled  by  their 
harmonic  number. 

Figure  9.  TEmni  mode  start  oscillation  condition  for  a  the  ridged  rectangular  cavity  for  an 
electron  beam  with  72  keV  and  0,  —  0.28.  All  modes  with  m  =  0, 1,  2  and  n  =  1, 
2, 3, 4, 5  are  shown.  The  arrow  shows  the  TEoai  fith  harmonic  emission  resonance. 

Figure  10.  Variation  of  the  start  oscillation  condition  for  the  third  harmonic  electron  cy¬ 
clotron  maser  interaction  (solid  curves)  and  second  harmonic  peniotron  interac- 
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tion  (dashed  curves)  as  a  function  of  dfa  for  the  TEsn  mode.  Note  that  in  the 
notation  used  here  that  d/a  =  i2o/i2.  The  beam  axial  velocity  \a  P,  =  0.1.  The 
cylindrical  ridged  cavity  parameters  are  $o  =  x/2JV  and  N  =  6. 
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Gyro-Pbniotron  Embsion  in  a  Magnetron- Type  Cavity* 

P.  VlTELLO 

Science  applications  International  Corporation 
MCLEAN,  VA  22102  USA 

ABSTRACT 

We  present  here  a  new  linear  theory  formulation  of  the  average  energy  loss  via 
the  electron  cyclotron  maser  and  peniotron  interactions  in  a  slotted  magnetron-type 
cavity.  Use  is  made  of  the  linear  theory  results  to  study  start  oscillation  conditions  and 
mode  competition  for  the  high-harmonic  gyro-peniotron  oscillator.  For  a  magnetron- 
type  cavity,  the  gyro-peniotron  is  found  to  have  the  same  sensitivity  to  guiding  center 
spread  as  for  a  smooth  walled  cavity.  Mode  competition,  as  well  as  the  start  oscillation 
beam  power  can  be  reduced  however  in  a  slotted  cavity  geometry.  In  the  non-linear 
regime,  the  saturated  emission  efficiency  for  a  magnetron-type  cavity  gyro-peniotron  is 
predicted  to  be  as  large  as  has  been  calculated  for  a  smooth  walled  cavity  device. 


Introduction 


The  gyro-peniotron  oscillator  is  a  device  with  the  potential  for  extremely  high  sat¬ 
urated  efficiencies  for  RF  emission.  Theoretical  modelingl^~*)  suggests  that  efficiencies 
(perhaps  as  great  as  80-90%)  can  be  obtained  for  electron  beams  with  Anite  velocity 
and  guiding  center  spread.  Even  though  the  peniotron  interaction  may  lead  to  higher 
saturated  efficiencies  than  the  electron  cyclotron  interaction,  in  the  weak  field  linear 
limit  the  cyclotron  maser  interaction  is  generally  very  much  stronger.  The  relative 
weakness  of  the  peniotron  interaction  can  lead  to  mode  competition  and  interierence 
firom  cyclotron  maser  modes  with  lower  start  oscillation  beam  powers.  Use  of  a  slotted, 
magnetron-type  cavity  presents  a  means  of  both  lowering  the  gyro-peniotron  start  os¬ 
cillation  current,  and  at  the  same  time  reducing  mode  competition  We  give  here  a 
new  formulation  of  the  magnetron-type  cavity  linear  theory  averi^e  beam  energy  loss 
(e./.  Referrace  (lOj),  and  discuss  mode  competition,  start  oscillation  conditions,  and 
the  behavior  of  the  non-linear  saturated  emission  efficiency  in  such  a  system. 

We  will  consider  a  magnetron-type  cavity  with  N  slots  each  of  width  20o.  The 
cavity  inner  radius,  outer  radius,  and  axial  length  are  respectively  a  and  b,  and  L.  For 
such  a  system,  the  beam  averaged  change  in  energy  6^  for  an  axis  symmetric  beam  of 
guiding  center  fg  is  given,  in  the  linear  regime,  by 
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where  Kg  =  k^ra  =  k^rit  u  the  perpendicular  wavenumber,  ri  is  the  Larmor 
radius,  F  =  ni  +  ji\r,  m  is  an  integer  which  determines  the  mode,  A;||  =  lif/L,  u  is  the 
cavity  frequency,  Bo  is  the  RF  field  amplitude. 
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Dimensionless  units  with  a  as  the  unit  of  length,  a/e  as  the  unit  of  time,  and  mgC^I\t\a 
as  the  unit  of  electric  and  magnetic  fields  are  used.  Each  term  in  j  in  above  relation  for 
6^  is  exactly  times  the  value  found  for  the  change  in  energy  for  a  smooth  walled  cavity 
with  the  single  mode  Electron  cyclotron  maser  emission  at  the  sth  cyclotron 

harmonic  comes  from  the  terms  containing  ^{x)  and  dx^p(z)/dx.  Electron  cyclotron 
absorption  at  the  sth  harmonic  comes  from  the  terms  in  p(x).  Peniotron  emission  at 
the  (s  —  l)th  harmonic  and  absorption  at  the  (s  +  l)th  harmonic  come  respectively  from 
the  terms  in  and  g{x~).  For  fixed  j  and  s,  the  relative  strengths  of  each  of  these 
interactions  are  the  same  as  in  the  smooth  cavity  case. 

If  the  beam  is  exactly  centered  on  axis,  only  the  harmonics  s  =  F  contribute  to 
Electron  cyclotron  emission  occurs  then  at  the  Fth  harmonics,  while  peniotron 
emission  occurs  at  the  (F  -  l)th  harmonics.  For  N  >2  there  is  no  interference  betwetm 
the  two  interactions  of  a  fixed  m  mode.  The  electron  cyclotron  maser  interaction  and 
peniotron  absorption  may  interfere  with  peniotron  emission  if  there  is  a  finite  spread 
in  guiding  centers.  As  the  relative  strengths  of  the  interactions  are  the  same  as  in  the 
smooth  cavity  case,  the  limits  on  guiding  center  spread  are  also  the  same  (see  Reference 
[11]).  Even  with  sero  guiding  center  spread,  mode  competition  may  take  place  between 
peniotron  emission  from  the  desired  Operating  m  mode  and  electron  cyclotron  maser 
emission  from  different  m  modes.  Such  mode  competition  can  be  a  serious  problem  for 
the  high-harmonic  smooth  walled  gyro-peniotron.I^I  A  N  slotted  magnetron-type  cavity 
cavity  has  an  imposed  iVth  fold  symmetry,  however,  which  is  lacking  in  a  smooth  walled 
cavity.  Operating  at  the  high-harmonic  m  =  F  =  a  =  N/2  optimises  the  amplitude  ap, 
leading  to  a  very  strong  beam  to  cavity  coupling.  For  m  modes  that  do  not  ‘^t”  well 
into  the  cavity  synunetry,  ap  is  quite  small  and  the  beam  to  cavity  coupling  is  weaker. 
The  overall  effect  is  a  reduction  in  mode  competition  between  the  m  =  IV/2  mode  and 
other  m  modes.  Besides  this  reduction  in  mode  competition,  a  slotted  cavity  can  also 
have  the  desirable  effect  of  moving  the  RF  field  radial  peik  inward  towards  the  beam. 
This  produces  a  stronger  interaction,  decreasing  the  start  oscillation  current,  which  is 
relatively  high  for  the  smooth  walled  gyro-peniotron.^^^l 

In  the  non-linear  limit,  as  in  the  linear  regime,  neighboring  s  harmonics,  differ  by 
at  least  by  for  a  beam  with  small  guiding  center  spread  about  the  axis.  If  iV  is  large 
and  is  one  operates  at  the  mode  m  =  N/2,  only  the  s  =  N/2  harmonic  will  show  a 
strong  coupling  between  the  beam  and  the  cavity.  The  other  harmonics  are  weaker  due 
to  having  either  much  smaller  amplitudes  ap  or  due  to  the  beam  particles  experiencing 
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very  large  phaM  drifts  relative  to  the  RF  field  during  their  passage  through  the  cavity. 
Since  only  one  harmonic  would  therefore  be  involved  in  gjrro-peniotron  magnetron- 
type  cavity  interaction,  we  e3q>ect  the  same  extremely  high  saturation  efficiencies  to  be 
obtmnable  as  have  been  calculated  for  smooth  cavity  systems. 

CONCLUSION 

We  show  from  a  linear  theory  analysis  of  a  slotted  magnetron-type  cavity  oscillator 
that  for  each  sth  harmonic  of  the  cyclotron  frequency,  the  relative  relative  strength  of  the 
electron  cyclotron  maser  interaction  and  the  peniotron  interaction  terms  are  the  same  as 
for  a  smooth  walled  cavity.  This  leads  to  a  gyro-peniotron  with  a  magnetron-type  cavity 
having  the  same  sensitivity  to  guiding  center  spread  as  is  found  for  a  smooth  walled 
system.  The  optimized  saturated  emission  efficiency  for  the  slotted  cavity  is  predicted 
to  be  also  as  high  as  has  been  calculated  for  the  smooth  walled  gyro-peniotron.  A 
slotted  cavity  design  is  likely  to  be  preferable,  however,  as  the  imposed  symmetry  can 
lead  to  a  reduction  in  mode  competition,  and  to  a  decrease  in  the  start  oscillation  beam 
current. 
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Start  oscillation  conditions  are  considered  for  slotted  rectangular  gyrotron  cavities. 
The  energy  loss  of  the  electron  beam  to  the  cavity  RF  field  is  presented  in  a  form  where 
the  geometric  and  the  gyrotron  instability  terms  at  each  harmonic  have  been  separated. 
It  is  found  that  the  addition  of  slots  to  a  cavity  lowers  the  required  beam  power  for  the 
start  oscillation  condition  and  can  lead  to  a  decrease  of  the  start  oscillation  beam  power 
as  the  harmonic  number  is  increased.  Previous  gyrotron  cavity  designs  have  consistently 
shown  the  start  oscillation  be^un  power  increasing  with  harmonic  number.  These  slotted 
devices  thus  appear  potentially  attractive  as  a  means  of  achieving  millimeter  wave 
emission  for  use  in  RF  heating  and  plasma  diagnostics. 


PACS  numbers:  85.10.Ka 


Rtni—iftn  in  n  gyrotrott  occun  at  the  harmonics  of  the  electron  cyclotron  frequency. 

t 

Promising  applications  for  high  povrer,  high  frequency  gyrotrons  include  plasma  heat¬ 
ing,  advanced  accelerators,  and  spectroscopy.  As  the  frequency  increases,  however, 
^  either  large  magnetic  fields  must  be  used  or  operation  at  high  cyciotron  harmonics 

is  required.  Presently,  moat  gyrotron  development  has  focused  on  devices  using  the 
first  or  second  harmonic.  These  low  harmonics  necessitate  either  the  use  of  a  super- 
^  conducting  magnet  or  a  pulsed  magnetic  field  if  millimeter  and  sub-millimeter  wave 

frequencies  are  to  be  achieved. 

High  harmonic  gyrotron  emission  has  been  demonstrated  for  several  tube  designs. 
Using  smooth-walled,  cylindrical  cavities  with  an  axis  encircling  beam,  harmonics  as 
high  as  m  =  11  have  been  observed  for  a  TEmii  ‘Svhispering  gallery”  mode  gyrotron 
oscillator.^'*  High  energy,  large  Larmor  orbit  beams  are  required  for  “whispering 
gallery”  mode  gyrotrons  as  the  modes  are  concentrated  towards  the  wall  of  the  cavity. 
Also  with  an  axis  encircling  beam  with  high  energy  (2  MeV)  electrons,  strong  emission 
at  the  12th  harmonic  has  been  generated  using  a  slotted  cylindrical  tube.*  Further 
investigations  of  slotted  rectangular  and  ridged  rectangular  tubes^~*  and  of  slotted 
cylindrical  tubes has  shown  that  even  for  moderate  beaun  energies,  strong  high 
harmonic  fringe  fields  can  be  positioned  at  the  beam  orbits.  All  of  these  studies  show 
'  that  slotted  structures  greatly  increase  the  high  harmonic  interaction  in  a  gyrotron. 

While  high  harmonic  gyrotron  oscillators  and  gyro-TWT  amplifiers  have  been 
studied  for  cylindrical  slotted  tubes,  only  gyro-TWT  amplifiers  have  been  been  inves- 

i' 

tigated  in  any  detail  for  rectangular  slotted  systems.  In  this  letter,  we  study  high 
harmonic  emission  in  rectangular,  slotted  gyrotron  oscillators.  We  find  that  the  addi¬ 
tion  of  slots  to  a  rectangular  gyrotron  cavity  greatly  enhances  the  interaction  between 

( 

the  electron  beam  and  the  cavity  RF  field  at  high  harmonics.  This  enhancement  is  so 
great  that,  under  some  circumstances,  the  interaction  strength  can  actually  mereose 
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with  hanumic  number.  As  a  consequence,  these  devices  are  a  potentially  interacting 
source  of  mUlimeter  wave  radiation. 

The  slotted  rectangular  cavity  considered  is  shown  in  Figure  1.  We  will  use  N  to 
designate  the  number  of  slots  and  Ir«  the  axial  length  of  the  cavity.  The  inner  and 
outer  slot  y  dimensions  are  R  and  while  the  cavity  width  in  the  x  direction  is  given 
as  £«.  The  distance  2d  between  the  slots  equals  twice  the  distance  from  the  first  or 
last  slot  to  the  tube  wall,  while  the  slot  width  is  2w.  We  expect  that  qualitatively 
our  results  would  apply  to  other  slotted  or  ridged  rectangular  cavity  designs  which 
show  strong  fringe  fields,  such  as  that  used  for  the  peniotron.^^^  In  the  following  we 
make  use  of  dimensionless  units,  with  the  cavity  transverse  length  R  being  our  scaling 
parameter.  In  these  units,  length  is  measured  in  units  of  i2,  frequency  in  units  of  c/i2, 
and  the  background  magnetic  and  RF  fields  are  measured  in  units  of  m,c^l\e\R^  where 
m«  is  the  electron  rest  mass,  e  is  the  electron  charge,  and  c  is  the  speed  of  light.  The 
dimensionless  cyclotron  frequency  O/'y  is  equal  in  these  units  to  B0/7  where  is  the 
assumed  constant  background  axial  magnetic  field  and  7  is  the  Lorentz  factor. 

To  determine  the  RF  fields,  it  it  convenient  to  expand  the  RF  field  separately  in 
the  cavity  proper  (Region  I)  and  in  the  slots  (Region  II)  in  terms  of  the  local  regional 
eigenfunctions,  and  to  then  match  the  fields  across  the  slot  openings.  We  take  the 
relative  phase  dependence  of  each  of  the  N  slots  as  being  cos[7r(9  -h  l/2)m//V],  where 
q  =  0, . . . ,  AT  —  1  is  the  slot  number.  The  value  of  the  mode  number  m  determines 
the  overall  mode.  The  2x  mode  corresponds  to  m  =  0,  while  the  it  mode  is  given  by 
m  ^  N.  la  order  to  obtain  a  closed  analytic  form  for  the  field  solution,  an  infinite 
sum  over  F,  where  kg  =  TirlLg,  is  used  only  in  Region  I,  with  just  a  single  term, 
kg  =  0,  being  kept  in  Region  n.  This  approximation  is  described  in  Ref.  7.  Only  the 
non-negative  terms  which  satisfy  F  =  2N}  ±  m,  and  j  =  0, 1,2, . . .  are  allowed.  The 
cavity  frequency,  w  =  (k^  +  kg)  for  the  /th  axial  mode  (kg  =  lirfLg)  is  determined 
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from  the  cavity  mode  m  by  solving  the  dispersion  relation 
2Nw  ^  k± 


(1) 


where  =  (**  - 

To  solve  for  the  interaction  of  the  beam  with  the  TE  RF  cavity  field  and  a  constant 
axial  magnetic  field  we  have  used  the  linearized  relativistic  single  particle  equation  of 
motion,  dU/dz  =  —{'i/U,)  ^E-\-(U  x  B/Tr)j  in  the  guiding  center  coordinates  (xge, 
ygc),  where  (f  =  la  the  product  of  the  electron  velocity  0  and  the  Lorentz  factor. 
After  averaging  over  initial  time  tg  and  phase  angle  tpg  of  the  beam,  the  change  in  the 
beam  energy  after  passing  through  the  oscillator  cavity  {S'))t,  can  be  calculated  using 
the  standard  techniques  of  linear  theory.^*  The  details  will  be  presented  elsewhere. 
We  find  that  can  be  expressed  in  the  compact  form 

3 

Y^ArPr,  [Re{f,)  +  »1m(/.)I,  (2) 

r 

where  Ap  is  the  Tth  amplitude  of  the  RP  field  in  Region  I,  Dp,  is  the  s  harmonic 
contribution  to  the  guiding  center  expansion  of  the  Fth  RF  field  term,  and  T,  gives 
the  s  harmonic  gyrotron  interaction  terms.  Again,  only  the  non-negative  terms  which 
satisfy  F  =  2Nj  ±  m,  where  j  =  0,1,2,...  are  allowed  in  equation  (2). 

The  full  behavior  of  is  given  by  equation  (2).  The  guiding  center  depen¬ 
dence  is  contained  entirely  in  the  factor  which  is  sensitive  to  the  fringe 

fields  generated  by  the  slots,  while  the  resonance  behavior  as  a  function  of  magnetic 
field  about  the  cyclotron  humonics  is  given  by  7,  which  contains  the  cyclotron  maser, 
Weibel,  and  peniotron  resonances.  The  contributions  from  Re(^)  and  Im(^)  re¬ 
spectively  determine  the  gyrotron  start  oscillation  beam  power  and  resonant  frequency 
detuning. 

Equation  (2)  can  be  simplified  if  a  ribbon  beam  parallel  to  the  x  axis  is  used. 
Averaging  over  the  x,e  guiding  center  position  leads  to  a  reduction  of  the  double  sum 


(^'y)uPo=  S 

— OO 
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over  cavity  eigenfunctions  inherent  in  |53r  =  2r  I3a  ^a^a»»  ^  single 

sum  if  the  ribbon  beam  is  uniformly  distributed  at  a  fixed  y^e*  After  averaging  the 
product  I]r -^f^rtllA^A^A*  replaced  by  a  factor  of  the  form 

Realistically,  however,  the  averaging  can  not  be  taken  from  Xffc  =  0  to  ^  as  the 
beam  cam  not  be  placed  closer  to  the  cavity  wall  than  its  Larmor  radius.  If  m  =  2Ni, 
where  t  =  0,1,...  then  there  is  a  AT-fold  periodicity  in  x  for  the  RF  field  and  the 
simplifying  guiding  center  averaging  can  equivalently  be  taken  over  auiy  of  the  regions 
from  Xge  =  Lx{j/N)  to  Xge  =  Lxik(N)t  where  0  <  j  <  iV  —  1,  and  j  <  k  <  N. 
This  also  implies  that  the  guiding  center  averaged  energy  loss  would  be 

unaffected  by  the  extension  of  the  cavity  in  the  x  direction  to  include  addition  slotted 
regions. 

The  cavity  geometry  determines  kx  and  the  magnitude  of  each  s  harmonic  con¬ 
tribution  to  the  overall  RF  cavity  field  of  a  particular  mode.  For  each  s  humonic, 
varies  with  the  cavity  geometry  only  through  changes  in  k±.  Since  the  addition  of  slots 
may  generate  strong  fringe  fields  without  a  significant  change  in  kx>  the  predominant 
change  at  high  harmonic  to  generally  comes  from  the  factor  |53r 

The  start  oscillation  beam  power  and  the  frequency  detuning  equation  for  steady 
state  operation  for  suiy  TE  mode  can  be  found  directly  from  Taking  the 

unperturbed  beam  energy  to  be  Pi,,  the  total  net  energy  transferred  from  the  beam  to 
the  cavity  fields  is  given  by  Pb{6'i)t„  v»«/('7o  “  Quite  generally,  the  starting  power 
can  be  given  by 


QPi,  =  -8.6  X  10®W(iro  -  l)//?e((^ kW,  (3) 


and  the  frequency  detuning  due  to  the  presence  of  the  beam  is 

_ 1  Im((g‘7)t 

u)o  ~  2QRe((fi7)t.»».)’ 
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where  W  is  the  tima  averaged,  volume  integrated  RF  field  energy  of  the  cavity,  and  Wo 
is  the  cold  cavity  frequency.  For  gyrotron  tubes,  the  quality  foctor  Q  is  due  mainly  to 
diffractive  losses  and,  with  an  accurate  estimate  of  its  value  obtained  from  cold  cavity 
tests,  provides  a  good  approjdmation  of  the  actual  start  oscillation  power  required.  By 
contrast,  the  linear  theory  values  for  the  frequency  detuning  often  differ  significantly 
from  the  detuning  under  normal  high  field  operating  conditions^’  and  should  used  with 
caution. 

We  find  that  the  addition  of  slots  can  not  only  enhance  enormously  the  high 
harmonic  interaction  and  hence  greatly  lower  QPh,  but  this  enhancement  at  fixed  beam 
position  can  inereaae  with  increasing  harmonic  number.  This  rapid  decrease  in  the 
start  oscillation  condition  with  increasing  harmonic  number  has  not  previously  been 
noted  for  any  gyrotron  device.  We  note  that  the  strongest  coupling  to  our  cavity  slots 
occurs  for  m  =  0  modes.  For  these  modes  the  following  parameter  range  was  found 
to  result  in  strong  high  harmonic  interaction  with  increases  with  harmonic  number: 
Isinfcxl”^  22  1  (implies  large  fringe  fields  generated  by  the  slots);  L,  ^  2N/n,  where  n 
is  the  transverse  y  mode  number  (all  kyj,  for  F  7^  0  imeiginaxy);  and  8in2irNwfLx  —  1 
(strong  beam  coupling  to  the  lowest  order  F  ^  0  fringe  field  term).  We  expect  proper 
beam  placement  in  other  cavity  designs  with  fringe  fields  similar  to  those  in  a  slotted 
rectangular  cavity  would  also  result  in  enhanced  high  harmonic  coupling  with  increasing 
harmonic  number. 

As  an  example  of  how  the  addition  of  slots  to  a  rectangular  cavity  gyrotron  oscil¬ 
lator  enhances  the  high  harmonic  emission  we  will  consider  the  start  oscillation  con¬ 
dition  QPh  for  the  cavity  and  beam  parameters  investigated  by  Han  and  Ferendeci^ 
for  a  slotted  rectangular  cavity  gyro-TWT.  For  the  cavity  mode  we  use  m  =  0,  /  =  1, 
and  the  second  y  transverse  mode  (i.e.  the  TEoai  mode),  which  Han  and  Ferendeci 
found  to  give  strong  sixth  harmonic  emission.  The  cavity  parameters  are  L*  =  2.186, 
L,  =  15,  Ho  =  1.198,  w  =  0.05466,  and  W  =  3.  For  the  beam  parameters  we 
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use  a  total  energy  of  70.82  keV,  and  0a  =  0.268,  which  correeponda  to  =  1.139, 
C/j.  3  0.45086,  and  Ua  —  0.30514.  The  electron  beam  is  centered  in  R^on  I  at 
^ge  =  ^s/2,  j/ge  =  Figure  2  shows  the  'variation  of  the  magnitude  the  RF 
electric  field,  |£|  =  {El  +  /  £«8in(k,s),  across  the  cavity  opening  for  these  pa¬ 

rameters.  Shown  in  the  figure  are  the  values  for  the  total  field  energy  per  unit  length, 
yHf  s  MtjLat  and  the  field  energy  per  unit  length  in  the  slots,  The  large  fringe 
field  observed  in  this  figure  couples  very  strongly  with  high  hannonic  eixuasion.  This 
coupling  will  be  significant  at  high  harmonics  since  with  rising  hannonic  number  s 
the  growth  of  the  the  Larmor  radius,  rr,  oi  s/9j./h>,  allows  the  beam  to  increasingly 
penetrate  the  fringe  fields.  In  Figure  3  we  show  QP^  as  a  function  of  magnetic  field  for 
multiple  harmonic  emission.  Harmonica  greater  than  the  sixth  are  not  shown  as  the 
beam  Larmor  radius  would  be  greater  that  one-half  the  cavity  width  R.  For  compari¬ 
son,  we  also  show  QP^  for  a  smooth-walled  rectangular  cavity.  For  the  slotted  cavity, 
harmonics  for  s  =  1-3  are  dominated  by  the  electron  cyclotron  maser  instability.  The 
fourth  harmonic  is  due  to  the  pcniotron  instability.  For  harmonics  s  >  5,  both  the 
electron  cyclotron  maser  and  peniotron  instabilities  strongly  contribute. 

As  the  sensitivity  to  guiding  center  is  of  great  importance,  we  have  given  in  Figures 
4a-4d  the  function  £*  =  IHr  “  3,4, 5,6.  Only  the  portion  of 

Region  I  which  is  at  least  one  Larmor  radius  from  the  cavity  walb  b  considered.  From 
equation  (3)  it  is  clear  that  QPj  a  (f  •)“’•  One  can  see  that  the  interaction  at  each 
harmonic  rapidly  increases  near  the  slot  openings  due  to  the  fringe  fields,  and  that  for 
s  =  3,4,5  there  are  values  of  y^e  for  which  there  b  little  dependence  on  Xge  of  the  beam, 
cavity-coupling.  These  values  would  be  ideal  for  the  placement  of  a  ribbon  beam  with 
an  extensive  z  width.  For  sixth  harmonic  operation  with  a  ribbon  beam,  the  cavity 
design  would  need  to  be  modified  to  reduce  the  rapid  variation  in  QP^  with  Xgc  in  the 
accessible  region  of  the  cavity. 
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For  our  examplo,  tho  strong  high  hsnnonk  coupling  which  increaass  with  hannAnic 
number  wu  observsd  even  when  the  beam  quality  was  decreased  by  ruing  a  guiding 
center  spread  Sxgc  ~  0.2,  Sy^c  ~  0.1,  and  an  axial  velocity  spread  60,  ^0.20,.  A 
Gaussian  axial  velocity  distribution  was  assumed  with  constant  beam  energy.  The  net 
affect  was  actually  to  decrease  QI^  as  the  interaction  strength  in  the  linear  regime  is 
a  rapidly  increasing  function  of  both  yg,  and  0~^.  In  general,  the  magnitude  of  the 
beam  quality  affects  will  depend  upon  the  details  of  the  gyrotron  dawign, 

In  conclusion  we  find  that  the  high  harmonic  interaction  observed  here  for  the 
slotted  rectangular  cavity  shows  great  potential  for  the  practical  development  of  a  high 
frequency,  low  magnetic  field  gyrotron  oscillator.  Very  high  power  sources  are  possible 
in  a  wide  cavity  with  multiple  slots  if  a  ribbon  electron  beam  is  used  to  reduce  space 
charge  limitations.  As  the  basic  gyrotron  instabilities  do  not  depend  upon  the  cavity 
geometry,  high  efficiency  emission  is  expected  for  the  slotted  rectangular  gyrotron,  es¬ 
pecially  if  the  system  is  optimized  for  emission  via  the  peniotron  interaction.*®  Further 
research  into  the  questions  of  mode  competition  and  non-linear  efficiency  needs  to  be 
carried  out  in  order  to  determine  how  besH  to  take  advantage  of  the  slotted  cavity  fringe 
field  interaction  for  gyrotron  devices. 

This  work  was  supported  by  the  AFOSR  under  contract  F49620-86-C-065. 
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Fignr*  Captions 


FIG.  1.  Cro68>section~of  the  slotted  rectsngolar  gyrotron. 

FIG.  2.  RF  electric  field  magnitude  across  the  slotted  rectangular  cavity  for  the  TE021 
mode.  Contours  are  shown  at  intervals  of  1/15  of  the  value  of  |£|. 

FIG.  3.  TEosi  mode  start  oscillation  condition  for  the  rectangular  canty.  The  solid  curves 
gives  QPi  for  the  slotted  cavity,  while  the  dashed  curve  gives  QPi  for  the  case 
without  slots.  Each  emission  resonance  is  labeled  by  its  harmonic  number  s. 

FIG.  4(a).  Surface  and  contour  plot  for  €*  for  a  =  3.  Contours  are  shown  at  intervals  of  1/15 
of  the  fn^viTtwitn  value  of  £*. 

FIG.  4(b)  Surface  and  contour  plot  for  £*  for  a  =  4.  Contours  are  shown  at  intervals  of  1/15 
of  the  maximum  value  of  £*. 

FIG.  4(c).  Surface  and  contour  plot  for  for  a  =  5.  Contours  are  shown  at  intervals  of  1/15 
of  the  maximum  value  of  £*. 

FIG.  4(d).  Surface  and  contour  plot  for  £*  for  a  =*  6.  Contours  are  shown  at  mterve.ls  of  1/15 
of  the  maximum  value  of  £*. 


FIG.  1.  Cross-section  of  the  slotted  rectangular  gyrotroh. 
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FIG.  4(b)  Surface  and  contour  plot  for  £*  for  a  =  4.  Contours  are  shown  at  intervals  of  1/15 
of  the  maximum  value  of  f  *. 


I 


Theory  of  High  Harmonic  Gyrotron  oscillators 
WITH  Slotted  cross  Section  Structure* 


p.  Vitello  and  C.  menyuk^ 


I 

»  , 

'  Science  applications  International  Corporation 

1710  GOODRIDGE  DRIVE 
McLean,  Virginia  22102  USA 

t 

I 


I 


I 


I 


*  This  work  was  supported  by  the  AFOSR  under  contract  F49620-86-C-065,  and 
by  the  University  of  California  at  Los  Angeles  under  Contract  400093  through  the 
I  U.S.  Army  Research  Office  Contract  ARO-DAAG  29-82-K-0004. 

^  Permanent  Address:  Department  of  Electrical  Engineering 

University  of  Maryland 
Catonsville,  MD  21228 
and 

*  College  Park,  MD  20742 


Submitted  to  the  IEEE  Special  Issue  on 
High  Power  Microwave  Generation 

49 


f 

ABSTRACT 

t  A  linear  theory  analysis  of  gyrotron  oscillators  with  slotted  cross  section  is  used 

to  calculate  the  net  change  in  beam  energy  In  our  new  formalism,  geometric 

factors  are  cleau'ly  distinguished  from  the  geometry-independent  harmonic  resonance 
I  terms  which  are  due  to  the  fundamental  electron  cyclotron  maser  and  peniotron  in¬ 

teractions.  This  separation  of  the  interaction  terms  from  the  geometric  factors  greatly 
simplifies  the  physical  analysis,  and  leads  to  a  very  compact  form  for  The 

I  theory  is  applied  to  slotted  rectangular  oscillators  (which  have  not  previously  been 

treated)  and  to  slotted  cylindrical  oscillators  to  show  that  a  unified  expression  can 
be  obtained  for  the  start  oscillation  condition.  In  sample  applications  of  our  theory, 
I  it  is  demonstrated  that  slots  lower  the  start  oscillation  condition  in  both  cylindrical 

and  rectangular  geometries,  and  can  lead  to  a  decrease  in  this  condition  as  harmonic 
number  is  increased  in  the  rectangular  geometry.  The  use  of  these  slotted  devices  thus 
I  appears  quite  attractive  in  the  millimeter  wave  regime.  We  also  find  that  the  peniotron 

interaction,  which  is  easily  identified  in  our  formalism,  may  be  very  strong  in  slotted 
cavities. 

I 

I 

I 
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I.  Introduction 

The  gyrotron  is  well  known  as  a  coherent  microwave  source  capable  of  generating 
unprecedented  power  levels.  Gyrotron  emission  occures  at  harmonics  of  the  electron  cy¬ 
clotron  frequency.  As  the  emission  frequency  increases  either  luge  magnetic  fields  must 
be  used  or  operation  at  high  cyclotron  harmonics  is  required.  Presently,  most  gyrotron 
development  has  focused  on  devices  using  either  the  first  or  second  harmonics  [1-3]. 
These  low  harmonics  necessitate  the  use  of  a  superconducting  magnet  if  frequencies  as 
high  as  the  millimeter  and  sub-millimeter  wave  regimes  are  to  be  achieved. 

High  harmonic  gyrotron  emission  has  been  achieved  for  several  tube  designs.  Using 
smooth  walled  cylindrical  cavities  and  an  axis  encircling  beam,  harmonics  as  high  as 
m  =  11  have  been  observed  for  a  TEmU  “whispering  gallery”  mode  gyrotron  oscillator 
[4-5].  A  high  energy,  large  Larmor  orbit  beam  is  required  for  “whispering  gallery” 
mode  gyrotrons  as  these  modes  tend  to  concentrate  the  RF  field  towards  the  wall  of 
the  cavity.  Also  using  an  axis  encircling  beam  with  high  energy  (2  MeV)  electrons, 
Destler  et  al.  [6]  have  reported  the  generation  of  a  strong  burst  of  microwave  radiation 
at  the  12th  harmonic  from  a  slotted  cylindrical  tube.  This  new  interactive  slotted 
structure  leads  to  an  excellent  high  harmonic  interaction  and  good  mode  selection. 
Further  investigations  of  slotted  cylindrical  tubes  [7-15]  and  of  slotted  rectangular  and 
ridged  rectangular  tubes  [16-19]  has  shown  that  even  for  moderate  beam  energies, 
strong  high  harmonic  fringe  fields  can  be  positioned  at  the  beam  orbits.  Two  slotted 
tube  geometries,  shown  in  Figure  1,  have  been  investigated  in  detail;  these  are  the 
cylindrical  (or  magnetron  type)  geometry,  and  rectangular  geometries,  both  slotted  and 
ridged.  As  there  is  no  real  distinction  between  slotted  and  ridged  tubes,  we  will  refer 
to  them  all  as  slotted.  Both  gyro-TWT  amplifiers  [7-10]  smd  gyrotron  oscillators  [ll- 
15]  have  been  studied  theoretically  in  slotted,  cylindrical  geometries.  For  rectangular 
geometries  only  gyro-TWT  amplifiers  have  been  studied  in  detail  for  slotted  cavities 


[16-19].  All  of  these  studies  show  that  slotted  structures  greatly  increase  the  high 
harmonic  interaction  in  a  gyrotron  relative  to  a  smooth  wall  design. 

In  this  paper  we  extend  the  modeling  of  slotted  gyrotron  oscillators  by  presenting 
a  linear  kinetic  theory  applicable  to  gyrotrons  of  both  rectangular  and  cylindrical  cross 
section  design.  In  our  analysis  we  calculate  the  net  energy  loss  of  the  beam,  (^7)t„ 
which  is  needed  to  determine  the  cavity  start  oscillation  condition  and  the  frequency 
detuning.  We  give  {6'l)toPa  ^  ^  harmonics  in  which  each  harmonic  contri¬ 

bution  consists  of  a  geometry-independent  factor  due  to  the  fundamentid  interactions 
and  a  geometry-dependent  factor  which  varies  from  device  to  device.  The  geometry- 
independent  factor  consists  of  clearly  distinguished  terms  due  to  the  electron  cyclotron 
maser  and  peniotron  interactions. 

This  paper  is  organized  as  follows:  In  §11,  we  derive  the  TE  mode  fields  in  either 
rectangular  or  cylindrical  slotted  cavities.  Two  specific  examples  are  presented  which 
are  then  used  in  the  remainder  of  the  paper.  In  §III,  we  present  our  linear  theory 
and  determine  the  start  oscillation  condition  and  frequency  detuning  of  the  oscillator 
cavity  modes.  In  §rV,  we  present  sample  calculations  which  shows  both  the  strong 
high  harmonic  content  of  slotted  cavities,  and  the  importance  of  peniotron  emission  for 
these  systems. 


n.  Properties  of  Slotted  Cavity  Fields 

Let  us  consider  the  behavior  of  the  TE  modes  in  a  tube  with  a  slotted  cross  section. 
If  the  axial  RF  xnagnetic  field  is  written  in  the  form 

B,  =  At(x,  y)/(«)e“*"*  (1) 

(where  f{z)  =  «**•*  for  an  amplifier  and  f(z)  =  sin(A;«^)  for  an  oscillator),  then  from 
Maxwell’s  equations  [20]  the  solution  for  At  will  be  an  eigenfunction  of  the  equation 

(Vr  •  Vr  +  =  0,  (2) 

subject  to  the  boundary  conditions  of  an  assumed  perfectly  conducting  tube  wall.  In 
(2).  *i  =  is  the  transverse  wavenumber  eigenvalue,  and  is  the  transverse 

Laplacian.  The  transverse  RF  electric  and  magnetic  fields,  Et  and  Bt  may  be  obtained 
from  Bg  as  follows, 

{z  X  VtB,)  ,  (3) 

Two  sets  of  orthogonal  eigenfunction  basis  sets  are  commonly  used  to  describe  At 
in  slotted  cavities.  If  the  cross  section  is  naturally  specified  in  a  Cartesian  coordinate 
system,  then  the  solution 

OO 

AT{x,y)  =  ^  Ar  +  are'**''*)  +  6re**»r*') ,  (5) 

r=-oo 

where  +  k^^  =  k^,  should  be  used.  For  a  cross  section  with  cylindrical  symmetry, 

AT(rJ)=  f;  >lrCr(*xr)e*^,  (6) 

r=-oo 

should  be  used.  In  (6),  Cr(fcx»')  =  +  or^r(^x’‘)i  where  Jp  “^d  ip  ^^re  Bessel 

functions  of  the  first  wd  second  kind.  The  constants  Ap,  <*ri  and  bp  depend  on  the 
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tube  geometry.  For  both  of  these  eigenfunction  basis  sets  the  coordinate  dependence 
is  separable.  The  total  time  averaged,  volume  integrated  RF  field  energy  stored  in  the 
fields,  “W,  can  in  all  cases  be  written  as 

W  =  ^  j  dV  (i-M* +B-B*y  (7) 

In  calculating  W,  any  static  magnetic  field  energy  is  not  included. 

We  now  consider  the  RF  cavity  fields  for  two  specific  cases,  a  slotted  co-axial 
cylindrical  cavity  as  shown  in  Figure  la  and  a  slotted  rectangular  cavity  as  shown  by 
Figure  lb.  We  will  use  N  to  designate  the  number  of  slots  and  Lm  as  the  axial  length  of 
the  cavity.  In  the  cylindrical  co-axial  cavity,  the  central  co-axial  radius  is  given  by  i2j, 
while  the  inner  and  outer  slot  radii  are  respectively  R  wd  Ro.  The  angular  width  of  the 
slots  is  taken  to  be  20<,.  In  the  rectangular  cavity,  the  inner  and  outer  slot  dimensions 
are  R  and  Rot  while  the  cavity  width  in  the  z  direction  is  given  as  Lg.  The  distance 
2d  between  the  slots  equals  twice  the  distance  from  the  first  or  last  slot  to  the  tube 
wallj  the  slot  width  is  taken  to  be  2w.  In  the  following  we  make  use  of  dimensionless 
units,  with  the  cavity  transverse  length  R  being  our  scaling  parameter.  In  these  units, 
length  is  measured  in  units  of  iZ,  time  in  units  of  R/c,  frequency  in  units  of  c/R,  and 
the  background  magnetic  and  RF  fields  are  measured  in  units  of  l\e\R,  where  mg 
is  the  electron  rest  mass,  e  is  the  electron  charge,  and  c  is  the  speed  of  light.  The 
dimensionless  cyclotron  frequency  0/^  is  equal  in  these  units  to  Bo/'T  where  Bo  is  the 
assumed  constant  background  magnetic  field  and  7  is  the  relativistic  factor. 

To  determine  the  RF  fields,  it  it  convenient  to  treat  the  cavity  proper  (Region 
I)  and  the  slots  (Region  11)  as  separate  expansions  which  are  matched  across  the  slot 
openings  [21].  We  take  the  relative  phase  dependence  of  each  of  the  N  slots  as  being 
g*2wmqfN  cylindrical  slotted  case  and  cos[7r(q  +  l/2)m/iV|  for  the  rectangular 

slotted  case,  where  j  =  0, . . . ,  iV  - 1  is  the  slot  number.  The  value  of  the  mode  number 
m  determines  the  overall  mode.  The  2ir  mode  corresponds  to  m  =  0,  while  the  x  mode 
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is  given  by  m  =  N/2  for  a  cylindrical  slotted  cavity,  and  by  m  =  iV  for  a  rectangular 
slotted  cavity.  In  order  to  obtain  a  closed  analytic  form  for  the  field  solution,  an  infinite 
sum  is  used  only  in  Region  I,  with  just  the  F  =  0  term  being  kept  in  Region  II.  The  RF 
fields  for  the  slotted  co-axial  cylindrical  cavity  are  then  given  approximateiy  by  [12] 


Region  I 


< 


f 


C 


Er  =  -  Eo^  E 

jss  — OO 


rCf(fcxr) 

kxr 


sin()k,^)c*(^®-'*'‘^ 


Ee  =  -iEo~  ^  A^Cr{kxr)sm(k,z)e'^'^^ 

''-L  y=-oo 

E,=  0, 


Br  = 

Be  = 

B,= 


Eo^  xf  cf  ^(fcir)  cos(/t,«)e’(^® 


j  =  -oo 
oo 


j=-oo 


rcf(fexr) 

*xr 


cos(fc,z)e*^*'® 
Bo  Y.  >«fCf(*xr)sin(fc,s)e<(^®-"*J, 


j=-oo 


(8) 

(9) 

(10) 

(11) 

(12) 

(13) 


f 
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Region  II 


Er  =  0, 

Ee  = 

Em  =  0, 

Br  =  coe(fc.z)e-"*, 

=  0, 

B.  =  £;„«‘^"*^/^X"C5^(jfe^r)8in(fc.a)e-‘‘*'*, 

where 

I  _  8m(rg„)  1 

r  Cf'(fcx)’ 

ij//  _  jl__L_ 

°  ^cr(ikx)’ 

and 

Cf  (fcxr)  =  Mkxr)  -  4(fc±^)yr(fci.r)/yAfcxi2i), 

C^"(*xr)  =  Mkj_r)  -  j;,{kj.Ro)Yo{kj.r)/Y’{kM- 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 

(22) 

(23) 


A  prime  represents  differentiation  with  respect  to  the  argiunent.  Due  to  the  cavity  sym¬ 
metry,  the  only  harmonics  which  contribute  are  F  =  m  +  jN,  where  j  =  0,  ±1,  ±2, . . . 
For  these  fields,  the  dispersion  relation  determining  ibx*  l^he  total  field  energy  W 


are 


N$o  ^  /smrtf<,y 

»r  V  ) 

;ss~O0  '  ' 


Cf(k±)  Ci^(k±) 
Cf'(k±)  Co"'(*x)’ 


(24) 
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and 

+  [Cr'(*x)]’  +  ^Cf'(tj.)Cf(kj.)  (l  +  ^)] 

^  [Cj'(k^H.)l’  -  [CJ'(*J]’ 

-  [Ci''(*x)]  ’  -  ic|''(*x)Cj'(tx)  (l  +  ^)  ]  |-  (25) 

The  slight  difference  between  W  as  given  in  (25)  and  the  form  presented  by  Chu  and 
Dialetis  [11-12]  stems  from  our  use  of  the  field  energy  density  E  •  E*  +  B  ■  B*  while 
Chu  and  Dialetis  used  2E’E*.  If  the  RF  fields  were  calculated  exactly,  with  an  infinite 
series  expansion  for  both  Regions  I  and  11,  then  the  energy  in  the  magnetic  fields  would 
equal  that  in  the  electric  fields  and  both  expressions  for  IV  would  be  identical. 

The  rectangular  slotted  fields  (see,  e.g.  [18]  for  a  derivation)  for  our  tube  are 

Region  I 


Ex  =  -t-^Ep  fcyr>lf  cos[kx^x)  sin(fcyry)  8in(fc,z)e  (26) 

^v-  *-^^o^kxi.Al‘Sin[kxi.x)cos{kj,j.y)sm{kxz)e~"^*,  (27) 

*x  r 

Ex  =  0,  (28) 

Bx  =  sinikx^x)  coe{ky^y)  cos(/f,a)e~*"*,  (29) 

''-L  r 

By  =  J3*»r>fr‘:os(*»r®)®>"(*yry)cos(A:,a)«~’"*,  (30) 

*-L  r 

Ba=  Eq  AI>  cos(fc,r^)  C08(fcyr  v)  sin(fes‘g)e~*'*'*»  (31) 

r 


57 


Region  n 


Ex  =  ^i-^Eo  cos{ir{q  +  l/2)m/N)A" 8in(A:xy)  sin(fc,*)e  *"*, 
Ey  =  0, 

E,^  0, 

Bx  —  0, 

By  =  —■^EoCos(jr(q  -h  l/2)m/N)AQ^  Bin(k±y)  coa[kgz)e~^'^* , 
Bn  =  EoCos{K{q  +  l/2)m/iV)^o^cos(A;xy)sin(fc,z)e“*"*, 

where,  fc*p  =  Tv  (Lx,  ky^  =  {k\  -  and 


4  = 


A"  = 


2A;j^  8in(feapti;)(-l)J' 

(1  +  5ro)  ^*r  ^vr  siniky^y 
kj.Lx 

Nsm{kxil-Ro)y 


Only  r  =  ±m  +  2N],  j  =  0,1,2,...,  with  F  >  0  we  to  be  included  in  the  sums.  For 
m  =  0,  each  term  for  F  =  2Nj  must  be  counted  twice.  For  our  model  of  a  rectangular 
slotted  cavity,  the  RF  field  in  Region  I  is  not  affected  by  the  presence  of  the  slots  for 
the  modes  m  =  Ni,  where  i  =  1,3,5,...,  and  has  the  same  form  as  in  a  simple  smooth 
walled  cavity.  For  all  other  rectangular  cavity  modes,  the  dispersion  relation  for  fcx 
and  volume  integrated  energy  W  can  be  written  as 
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m.  Linear  Kinetic  Theory 

To  solve  for  the  interaction  of  the  electron  beam  with  the  RF  field  we  use  the 
single  particle  relativistic  equation  of  motion 


dz 


U, 


f? 

E+ - 


(42) 


where  U  =  'fj3  ]a  the  product  of  the  electron  velocity  and  the  Lorentz  factor  7.  Note 
that  we  are  using  time,  <(2),  as  a  dependent  variable,  and  the  axial  position  z  as  the 
independent  variable.  The  use  of  z  facilitates  the  comparison  of  the  linear  theory  with 
the  weak  field  limit  of  a  nonlinear  numerical  modeling  of  the  equations  of  motion.  In 
dimensionless  units,  ^  =  {1  +  U  •  0)^^^  =  {1  +  U,  +  gives  the  electron  energy. 

From  (42)  it  follows  that  the  chauige  in  'y  is 


d'^  _  U’E 
dz~  U,  ' 


(43) 


Regardless  of  the  symmetry  of  the  overall  cavity  fields,  the  guiding  center  coor¬ 
dinates  (see  Figure  2)  are  the  natural  system  to  use  when  carrying  our  linear  theory. 
The  guiding  center  variables  may  be  expressed  as  Ug,  <p,  Xgc,  y^o  0  Cartesian 
coordinates  or  (C/'±,  17,,  v?,  »‘jc>  V’yo  0  in  cylindrical  coordinates.  The  guiding  center 
coordinates  are  related  to  cavity  frame  coordinates  by  the  relations 


x  =  Xgc  +  ri,  cos(v?) 

=  Tgc  cos(V»pc)  +  rt  cos((p) ,  (44) 


I 


I 


> 


V  =  t/gc  +  rt  sin(^) 

=  rgcsm{rf>gc)  +  ri,  sin(v3). 

(45) 

If*  =  -t;j.sin(v>), 

(46) 

Uy  =  U±cos{ip). 

(47) 
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The  position  (xj*,  y^e)  or  (rye,  V'jc)  corresponds  to  the  electron  guiding  center  position, 
with  =  U±/Bo  being  the  Larmor  radius.  The  dynamical  equations  for  the  guiding 
center  variables  foHow  from  (42,  44-47), 


1 

iz-  + 

(48) 

iz  =  +  r.' 

(49) 

1 

dU.  Uj. 
dz  ~ 

(50) 

II 

Xi  1 

(51) 

1 

and  either 

^  +  V’B,,  h)  cos(v) 

t 

+  (Ert,  -  UgB^/'i  +  UxB,li)  sin(¥>)] , 

(52) 

1 

or 

-  (£?rt  -  U,B^h  +  UxB,li)  cos(ip)] , 

(53) 

1 

dz  ^  U,Bo  cos(0yc  -  'p) 

-  (Ert,  -  UgB^I'i  +  UxB^li)  sin(0yc  -  <p)\ , 

(54) 

1 

+  (^rt  -  UgB^f'i  +  C^x5,/7)cos(0yc  -  >p)] . 

(55) 

1 

The  RF  fields  terms,  E^,  Er^ ,  and  Br^,  correspond  to  the  values  of  S  and  B  along 

the  ^  and  fi  directions. 
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To  solve  for  the  perturbations  induced  by  the  RF  field,  we  linearize  the  equations 
of  motion  as  follows 


t 

U±  =  Uu  +  SU^, 

(56) 

u,=-u,^  +  su„ 

(57) 

1 

<P  =  <Po  +  Bo^  +6<p, 

(58) 

Xgc  —  Xgc„  +  SXge, 

(59) 

1 

Vge  —  VgCo  +  ^Vget 

(60) 

fgc  —  ^gco  "h  ^^gci 

(61) 

- 

^ge  —  ^gca  +  S4>gei 

(62) 

,  =  + 

(63) 

where  the  constants  <Poy  Xgco,  represent  the  initial 

variable  values  at  the  entrance  position  of  the  cavity  Zg.  In  the  following  we  take  Zo  =  0. 
The  equations  of  motion  for  the  perturbed  variables  then  follow  as 

(64) 

(65) 

(66) 

(67) 


-f- 

I' 

i 


dtft/x 

—  —  E° 

dz 

ddyj 

-  E 

dz 

ddC^, 

dz 

_  ^j-o  pQ 

-  vZ 

d^t 

_  1  i 

dz 

e 

o 

t- 

1 

1 

o 

tl 


6U, 
Ui  ’ 


1  +  ^L 
u.. 


6U.  -  U±JUu 
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X 


I 


I 


I 


I 


I 


I 


dSXgc 

-  ^0  \ 

dz 

U,„Bo  L 

dSygc 

-  \ 

dz 

-  U,^B,  L 

d^>’ac 

-  \ 

dz 

u.„B,  L 

dS^gc 

To 

dz 

U.^Bo 

gco 


where 


(68) 

(69) 

(70) 

(71) 


=  +  (72) 

Et,  =  Et,  -  U^^B%ho  +  VlB%ho.  (73) 

<l>o  =  <fio  +  Boz/UM,.  (74) 


E^,  E°^,  B^,  B°^,  arm  B°  are  the  RF  field  terms  in  the  tp,  fi,  and  I  directions 
evaluated  along  the  unperturbed  orbits.  E^  and  E°^  represent  respectively  the  <p  and 
fi  components  of  vector  E-\-U  'X.B l~i  from  the  equation  of  motion  also  evaluated  along 
the  unperturbed  orbits.  The  variable  <^o  gives  the  phase  variation  about  the  guiding 
center  for  zero  RF  field. 

To  calculate  the  average  change  in  energy  for  the  beam  we  must  also  linearize  the 
energy  equation  (43),  which  becomes 


d57 

dz 


6U^ 


?0  * 


+ 


VI 


’E^*- 


Ux. 


*o 


+  Sr 


dt 


gc 


dr, 


gc 


gc 


dtj), 


^  dE^*  ^  dE^^ 
+  Sri,-^  ^  +  Sip — ^ 


gc 


drt 


dp  J 


(75) 


or 


dS'^ 

dz 


■ 


O.  ,  UxJU.^o.  Ux„ 

T  ,rO 


U,„ 


dE^* 

St  — +  Sx, 


dt 


dE^ 

dXgc 


^  dE^/  ^  dE^*  ^  dE^*-' 
+  ^yge  +  ^'P - 


dygc 


drt 


dp  J 


(76) 
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For  a  constant  background  magnetic  field  Bo,  Sri,  =  SUxjBo-  The  equations  for 
dj7/dz  have  been  expanded  to  second  order  in  the  RF  field  amplitude  Eg. 

In  order  to  solve  these  linearized  equations  we  must  express  the  RF  fields  in  terms 
of  the  guiding  center  coordinates.  For  the  cylindrical  geometry  eigenfunction  set  (8-13) 
Graf’s  theorem  [22], 


-  0)  =  ^  Ci,+ki<gc)Jk{KL)^^k{ir -il}gc  +  <p),  (77) 


A:=— oo 


for  Tgc  >  ri,oT 


Cu{kxr)l°^t/{tl}gc  -  6)  =  Y",  '^^+k{<gc)Ck{<L)l^k[ir -^gc  +  ip),  (78) 


fc=-CX5 


for  fge  <  ri,  is  used  where  Kgc  =  k^ryc  and  Ki,  =  kx^L-  For  the  rectangular  geometry 
eigenfunction  set  (26-31),  the  expansions 


cos  [Ki  sin(v3  ±  A)]  =  ^  J*(Kl)  cos  [fc(^  ±  A)] , 

k=—oo 

oo 

sin  [Kisin(ip  ±  A)]  =  E  Jk{KL)sin  [k(<p  ±  A)], 

fc=:  — OO 
OO 

cos  [Ki  cos(y?  ±  A)]  =  Y2,  •7fc(Ki)  cos  [fc((p  ±  A  +  7r/2)], 

Jk=— OO 

oo 

sin  [Ki  cos(v9  ±  A)]  =  E  •7fc(Ki)  sin  [A:((p  ±  A  +  7r/2)] , 

fc=— oo 


are  similarly  used.  In  (79-82),  cos  A  =  k^/kx- 


Using  the  above  expansions,  the  guiding  center  RF  fields  in  Region  I  which  contains 


« 
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the  beam  can  be  expressed  as 


II 

(83) 

II 

r  *=-00 

(84) 

Br,= 

-^oE^r  E  i’r.J:(X:i)cos(k,z)e*(«»’-"*), 

r  «=-oo 

(85) 

(86) 

oo 

(87) 

r  «=— oo 


where  the  function  Drt  is 

=  (88) 

for  the  cylindrical  coordinates  eigenfunction  expansion,  and 

Dr,  =  J  *•«  cos  {ky^Vgc  -  «Ar)  +  cos  [ky^Vgc  +  «Ar)]  (89) 

for  the  Cartesian  coordinates  expansion.  For  cylindrical  coordinates,  if  ri,  >  Vgc,  the 
replacement  —*  J  and  J  must  be  made  in  (83-88).  The  form  of  the  guiding 
center  RF  fields  is  clearly  independent  of  the  choice  of  the  coordinate  eigenfunction 
basis  set  one  wishes  to  use  to  describe  the  cavity  fields  and  follows  from  the  Fourier 
transformation  for  our  separable  eigenfunction  expansion.  To  determine  the  contribu¬ 
tion  of  the  cavity  field  at  the  sth  harmonic,  one  merely  sums  the  contribution  from 
each  r  eigenfunction  AlPr,- 

When  evaluated  along  the  unperturbed  particle  orbits  the  only  variation  in  the 
guiding  center  RF  fields  (83-87)  comes  from  sin(fc*z)e*^*^“"*),  which  takes  the  sim¬ 
ple  form  («**'*  —  «“**'*•)  making  the  linearized  equa¬ 

tions  of  motion  simple  to  evaluate.  These  equations  are  solved  in  the  following  order: 
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We  first  integrate  (64,  66,  68-69)  or  (64,  66,  70-71)  to  determine  6U±_,  SUg,  and  either 
Sxge  and  Sygc  or  Stge  amd  The  expressions  for  6U±_  and  6Ug  are  then  substituted 
into  (65,  67)  to  determine  Stp  and  6t.  Using  these  perturbed  orbital  variables,  (75)  or 
(76)  is  finally  integrated  from  2  =  0to2  =  Lin  order  to  calculate  S').  The  perturbed 
energy  change  is  then  averaged  over  to  ^d  the  phase  angle  <po.  We  will  not  average 
over  guiding  center  position  at  this  point.  The  averaged  form  for  which  comes 

from  a  straight  forward  integration  of  (75)  or  (76)  is  given  by 


IC  \ 


r  A 


*=— oo 


(90) 


where 


tti  - 


(w*  -  ik^)  [t>_  (e*>-  +  1)  -  2  (c’>-  -  1)]  Hvaz 


{ifi-  -  e’**-  +  l)  ((w  -  kgUg^l'io)  Ttk,  -  kiUx„UTK»ho\  ,  (91) 


VI  Ll 


-  [;»•-«"-  -  +  >1  »r.. 

'loUg^unt  '■ 

-  — ^  («’■'*-  -  1)  ((W  +  k.U.J'io)  Ttk,  -  kxUuUvKzho] ,  (92) 

WM+M- 


Hr\* 


M±  —  (w  ~  sBof'io  i  kgUg^f'^o)  i'*'7p/^*o> 


(93) 
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and 


(94) 

Tvk.  =  Hvk.  +  <lJ:  (KlJ 

Ur  A.  =  (/CrJ  [IVViPa.-i  -  PM 

+  *^*+1  (^io)  [^r*+i^A*+i  ~  Pr»pA.$]  (96) 

The  functions  as  and  as  are  related  to  ai  and  04  by 

as  (*,)  =  ai  (-*,) ,  (97) 

as  (**)  =  04  (-fc,) .  (98) 

Equation  (90)  for  contains  a  double  sum  over  the  cavity  moae  harmonics  F 

and  A.  For  the  cylindrical  coordinates  expansion,  F  and  A  are  restricted  to:  F  =  Nj +m 
for  j  =  0,  ±1,  ±2,  . . .  and  A  =  Nk  +  m  for  fc  =  0,  ±1,  ±2,  . . ..  For  the  Cartesian 
coordinates  expansion,  the  non-negative  harmonics  which  satisfy  F  =  2Nj  ±  m,  where 
j  =  0, 1, 2, . . .,  and  A  =  2Nk  ±  m,  where  A  =  0, 1, 2, . . .,  are  allowed.  When  m  =  0,  the 
two  Cartesian  F  cavity  mode  harmonics  F  =  2Nj  +  m  and  F  =  2Nj  —  m,  and  the  two  A 
cavity  mode  harmonics  A  =  2Nk  m  and  A  =  2Nk  —  m  are  identical.  In  this  case,  the 
sums  over  F  or  A  can  be  written  as  twice  the  sum  for  all  F  =  2Nj  or  A  =  2Nk.  Note 
again  that  for  the  cylindrical  expansion,  if  ri,^  >  r^c.  one  must  make  the  replacement 
-*  J  and  J  — ♦  in  the  above  relations. 

In  (90),  the  variation  of  with  guiding  center  position  only  enters  through 

the  factors  ^r«^A«>  ^r«+i^A«+i»  and  pT»-\PKa-\‘  Averaging  over  the  guiding  cen¬ 
ter  position  V’bc,  from  0  to  x  for  the  cylindrical  case  or  Xgc„  from  0  to  L,  for  the 
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Cartesian  coordinate  case  collapses  the  double  sum  over  F  and  A  to  a  single  sum 
with  r  =  A  and  replaced  by  either  or  [co8(Aypyjc,)^co8(sAr)^  + 

sin(Aypy0c^)^sin(sAr)^]/2.  For  the  cylindrical  slotted  gyrotron  oscillator,  the  resulting 
expression  for  agrees  with  the  results  of  Chu  and  Dialetis  [11-12].  While 

this  guiding  center  averaging  is  reasonable  for  modeling  am  axis-symmetric  beam  in  a 
cylindricad  cavity,  it  u  not  reasonable  in  a  rectangulair  cavity,  where  Xgc^  camnot  re¬ 
alistically  be  taken  within  rt^  of  the  cavity  waJls.  If  m  =  2JVt  where  »  =  0, 1,2, ... , 
then  the  RF  field  shows  a  i\r-fold  periodicity  in  the  z  direction  and  the  simplifying 
Xgc„  averaging  can  be  done  from  Lj,{j/N)  to  Lt{kfN),  where  0  <  j  <  N  -  1,  and 
j  <  k  <  N. 

The  form  of  given  in  (90)  is  a  very  complicated  function  of  the  beam  and 

cavity  pairauneters.  It  is  far  from  the  most  compact  form  for  the  net  energy  change,  and 
it  does  not  clearly  show  the  separate  contributions  from  the  electron  cyclotron  maser 
interaction,  the  peniotron  interaction,  and  other  interactions.  A  more  compact  form 
can  be  obtained  by  combining  and  re-arranging  the  terms  in  01  +  02  +  03+04.  To 
start  off  with,  we  will  re-write  in  the  form 


)toV>o  ~ 


-L  r  A 


k±U. 


MiLm  kjiUz  \ 

4u;  lo  J  J 


where 


9  = 


g-t(/irA+(J-l)irl  ^  J 

2(A*-1)“ 


lirX 


(100) 

(101) 
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While  (99)  shows  explicitly  the  rapid  variations  in  the  magnetic  mistuning  variable 
A  about  the  harmonic  resonances,  one  can  still  obtain  a  simpler  form  for 
Upon  expanding  and  re-writing  the  JJrAo  ^rA«i  and  Uta*  terms,  the  following  final 
expression  can  be  obtained 


ih) 


toVo  ~  ~ 


oo  r 


\  'l„v^  V*.  u  dX  J 

- (<!,.)  »<*> 

\  '1o<d>  'lofd}  ) 

-  Kii  -  J.'L,  (<i.)  J(A-)] 

.XmL,  r j.”  (JCi.)  /  w  iy(A)  t.  <iAV(A)  t.  \ 

*  4  '1oUa„  \ka  dX  U)  dX  <JJ  ) 

-  (Kl.)  ^'^0 

-  [‘j.t'j..  Wli  («£..)  /(A-^)  -  /(A-)) 


+  k.u^  (J.+1  ^;a..  (>Ci.)  +  J..,  (IC,  )  J',  (>CiJ  +  tJ,  (K.)  j:  (Kl. 


where 


/  cos(|*)  \ 

’-V(A“-I)/  ’ 

(103) 

A  4sin(M)cos(/i) 

(A»-l)  '  *.i.(A»-l)>' 

(104) 

\  'lo  } 

(105) 

n  =  /ff  A/2  +  (/  -  l)x/2. 

(106) 

As  the  product  ^r#  I^A -^A^A* ‘®  real,  c..id  since  all  of  the  variables  and  functions 

within  the  brackets  {. . .}  are  real,  the  first  set  of  terms  in  (102)  give  the  real  part 
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of  the  second  set  of  terms  multiplied  by  i  give  the  imaginary  part  of 

As  is  shown  below,  and  Ini((d'y)t,^,)  respectively  determine 

the  gyrotron  start  oscillation  beam  power  and  resonant  frequency  detuning.  Averaging 
over  the  guiding  center  angle  tfige,  for  a  cylindrical  cavity  or  over  Xye,  for  a  rectangular 
cavity  again  leads  to  a  reduction  of  the  double  sum  over  cavity  eigenfunctions  to  a  single 
sum.  After  averaging  the  product  Er  Ea  “  replaced  by  Er 

where  lfr«P  is  either  (or  [J/.*]*  if  tl,  >  rjc,)  or  [cos(fcyj.yjc„)^  cos(sAr)^  + 

sin(kffp!/ge„)^ sm(sAr)^]/2  respectively  for  cylindrical  and  Cartesian  geometries.  The 
smooth  walled  limit  of  equation  (102)  is  obtained  by  taking  the  limit  w  0  or  »  0 
in  the  dispersion  relation,  and  by  using  F  =  m,  ^^^af  =  !>  =  0. 

The  full  behavior  of  is  clearly  given  by  equation  (102).  The  guiding  center 

dependence  is  given  by  the  factor  Af  the  resonance  behavior  as  a 

function  of  magnetic  field  about  the  cyclotron  humonics  is  given  by  the  functions  y 
and  /;  and  the  several  interactions  which  simultaneously  take  place  in  a  gyrotron  are 
explicitly  separated.  Cyclotron  maser  and  Weibel  emission  resonances  are  given  by 
the  tetms  in  •  •)•  The  terms  in  2sJ'(. ..)  are  due  to  cyclotron  maser 

absorption.  Peniotron  emission  and  absorption  resonances  are  represented  respectively 
by  the  terms  with  the  factors  ff(A~),  /(A“),  ff(A'''),  and  /(A'*').  We  note  that  equation 
(102)  for  a  smooth  walled  cylindrical  cavity  is  equivalent  to  the  relation  derived  by 
Brand  [23],  and  that  the  real  portion  of  {S'^)t„<po  averaged  over  is  in  the  same 
form  given  previously  by  Vitello  [24].  For  slotted  cylindrical  oscillators,  the  real  part  of 
equation  (102)  can  also  be  shown  to  equal  the  more  complicated  relation  by  Chu  and 
Dialetis  [11-12]. 

The  start  oscillation  beam  power  and  the  frequency  detuning  equation  for  steady 
state  operation  for  any  TE^ni  mode  can  be  found  directly  from  (^'y)tay>o'  Taking  the 
unperturbed  beam  energy  to  be  the  total  net  energy  transferred  from  the  beam 


to  the  cavity  fields  is  given  by  Ph{6'l)uv»l{'1o  —  !)•  Quite  generally  [25],  the  starting 
power  can  be  given  by 

QPt  =  -8.6  X  10« V  ('lo  -  1)/Re((«7)t.  ^ J  kW,  (107) 


and  the  frequency  detuning  due  to  the  presence  of  the  beam  is 

_ _ 1 

where  Wo  =  (k*  +  fcj)  is  the  cold  cavity  frequency.  For  gyrotron  tubes,  the  quality 
factor  Q  is  due  mainly  to  diffractive  losses,  and  with  a  good  estimate  of  its  value 
obtained  from  cold  cavity  tests  (107)  provides  an  good  estimate  of  the  actual  start 
oscillation  power  required.  By  contrast,  the  linear  theory  values  for  the  frequency 
detuning  often  differ  significantly  from  the  detuning  under  normal  high  field  operating 
conditions  [23]  and  should  used  with  caution.  Again  we  stress  that  for  rt,^  >  rge„  there 
must  be  a  replacement  in  the  above  relations  of  J  amd  J  -*  C^. 


70 


IV.  Discussion 


We  have  extended  the  modeling  of  the  linear  kinetic  theory  for  gyrotron  oscillators 
to  cavities  with  slotted  rectangular  and  cylindrical  co-axial  cross  sections.  The  averaged 
beam  perturbed  energy  change  is  shown  in  a  form  (102)  which  is  the  same  in  rectangular 
and  cylindrical  cavities.  The  cavity  geometry  determines  only  k±  and  the  magnitude 
of  each  harmonic  contribution  to  the  overall  RF  cavity  field  of  a  particular  mode.  For 
each  harmonic,  the  functional  form  of  the  differing  interactions,  such  as  the  electron 
cyclotron  maser  and  the  peniotron  interactions,  do  not  vary  with  the  cavity  geometry 
and  are  the  same  for  smooth  walled  and  slotted  tubes.  Studies  of  the  relative  strengths 
of  the  interactions  for  smooth  walled  gyrotrons  [25]  can  therefore  be  applied  directly 
to  slotted  cavity  devices. 

We  apply  here  our  results  to  the  slotted  rectangular  and  cylindrical  cavities  dis¬ 
cussed  earlier.  In  the  case  of  a  cylindrical  cavity  we  set  R,-  =  0,  average  <(>„  over 
000,  >  use  an  axis  encircling  electron  beam  to  allow  for  comparison  with  the  results 
of  Chu  and  Dialetis  [11-12].  In  Figure  3  we  show  the  £$  RF  field  profile  for  the  TEsu 
mode  for  a  cavity  with  iV  =  6  and  Ro  —  1.4.  As  one  of  the  major  advantages  of  the  slot¬ 
ted  cylindrical  cavity  is  the  reduction  of  the  start  oscillation  condition  with  increasing 
Ro,  we  plot  QP^  in  Figure  4  for  both  the  electron  cyclotron  maser  interaction  and  the 
peniotron  interaction  as  a  function  of  Ro.  In  the  figure,  Lg  is  increased  with  increasing 
Ro  to  hold  ifex/ik||  fixed  at  10.7.  The  decrease  in  QP\,  with  Ro  is  partially  due  to  the 
effects  of  a  longer  cavity,  but  is  primarily  caused  by  the  increasing  slot  depth.  The 
electron  cyclotron  maser  emission  curve  for  49.9  KeV  is  terminated  at  the  Ro  value 
where  the  Larmor  radius  strikes  the  cavity  wall.  Our  results  for  the  electron  cyclotron 
maser  QPi,  are  in  agreement  with  Chu  and  Dialetis.  However,  we  find  that  at  low  beam 
energy  the  start  oscillation  condition  is  actually  smaller  for  the  peniotron  interaction. 
Chu  and  Dialetis  give  results  only  for  the  electron  cyclotron  maser  resonances  amd  do 
not  consider  the  peniotron  resonances. 


As  an  example  of  the  effect  of  slots  in  a  rectangular  cavity  gyrotron  oscillator  we 
will  consider  the  start  oscillation  condition  QP^  for  the  cavity  and  beam  parameters 
investigated  by  Han  and  Ferendeci  [17]  for  a  slotted  rectangular  cavity  gyro-TWT.  For 
the  T£o3i  mode,  which  Han  and  Ferendeci  found  to  give  strong  sixth  harmonic  emis¬ 
sion,  we  place  the  electron  beam  in  Region  I  at  Zpe,  =  I'*/2,  ygc„  =  R/2.  Figure  5 
shows  the  Ex  component  of  the  RF  field.  The  strong  fringe  field  observed  in  this  figure 
couples  very  strongly  with  high  harmonic  emission.  This  coupling  will  be  significant 
at  high  harmonics  since  with  rising  harmonic  number  s  the  growth  of  the  the  Larmor 
radius,  ri,^  allows  the  beam  to  increasingly  penetrate  the  fringe  fields.  In 

Figure  6  we  show  QPb  for  the  slotted  rectangular  cavity  as  a  function  of  magnetic  field 
for  multiple  harmonic  emission.  Harmonics  greater  than  the  sixth  are  not  shown  as  the 
beam  Larmor  radius  would  be  greater  that  one-half  the  cavity  width  R.  For  comparison 
QP(,  for  a  smooth  walled  rectangular  cavity  is  shown  also.  For  our  beam  position,  the 
smooth  walled  cavity  only  shows  emission  at  the  even  harmonics.  For  the  slotted  cavity, 
harmonics  for  s  =1-3  are  dominated  by  the  electron  cyclotron  maser  instability.  The 
fourth  harmonic  is  due  primwily  to  the  peniotron  instability.  For  harmonics  greater 
than  fourth,  both  the  electron  cyclotron  maser  and  peniotron  instabilities  strongly  con¬ 
tribute.  It  is  evident  that  the  addition  of  slotted  not  only  enormously  enhances  the 
very  high  harmonic  (s  >  4)  interaction  and  hence  greatly  lower  QPb,  but  this  enhance¬ 
ment  increases  with  increasing  harmonic  number.  This  decrease  in  QPb  with  harmonic 
number  is  implicit  also  in  the  ridged  gyro-TWT  modeling  of  Hzm  and  Ferendeci.  For 
the  smooth  wailed  rectangular,  smooth  walled  cylindrical,  or  even  ridged  cylindrical 
cavities  with  axis  encircling  beam  previously  studied  in  the  literature,  the  stand2ird 
behavior  with  increasing  harmonic  number  is  a  growth  in  the  start  oscillation  beam 
power.  The  high  harmonic  interaction  observed  here  for  the  ridged  rectangular  cavity 
shows  a  great  potential  for  the  practical  development  of  a  high  frequency,  low  magnetic 
field  gyrotron  oscillator. 


In  conclusion  we  find  that  for  complex  cross  gyrotron  cavities  the  start  oscillation 
condition  can  be  written  in  a  form  which  separates  the  geometric  factors  from  the  phys¬ 
ical  interaction  terms.  An  analysis  of  slotted  cavities  using  this  formalism  shows  that 
slots  can  be  used  to  enormously  enhance  the  high  harmonic  RF  field,  beam  interaction 
for  both  rectangular  and  cylindrical  cavities,  and  for  both  the  fundamental  electron 
cyclotron  maser  and  the  peniotron  interactions. 


References 

[1]  J.M.  Baird,  "Survey  of  fast  wave  tube  developments,”  Tteh.  Dig.  Ini.  Electron 
Devices  Meet.,  pp.  156-163,  1979. 

[2]  R.S.  Symons  Sc  H.R.  Jory,  "Cyclotron  resonance  devices,”  in  Advances  in  Elec¬ 
tronics  &  Electron  Physics,  L.  Marton  Sc  C.  Marton,  Eds.,  vol.  55,  New  York: 
Academic,  pp.  1-75,  1981. 

[3]  V.L.  Granatstein,  M.  Read,  k  L.R.  Barnett,  "Measured  performance  of  gyrotron 
oscillators  and  amplifiers,”  Int.  J.  Infrared  Millimeter  Waves,  vol.  5,  pp.  267-304, 
1984. 

[4]  D.B.  McDermott,  N.C.  Luhmann,  Jr.,  Sc  D.S.  Furuno,  "Operation  of  a  high- 
harmonic  gyrotron,”  Eighth  Int.  Conf.  Infrared  Millimeter  Waves,  TH4.1,  1983. 

[5]  D.B.  McDermott  k  N.C.  Luhmann,  Jr.,  "Operation  of  a  compact  mm-wave  high- 
harmonic  gyrotron,”  Proc.  SPIE,  vol.  423,  pp.  58-63,  1983. 

[6]  W.W.  Destler,  R.L.  Weiler,  k  C.D.  Striffler,  “High  power  microwave  generation 
from  a  rotating  e-layer  in  a  magnetron-type  waveguide,”  Appl.  Phys.  Lett.,  vol.  38, 
pp.  570-572,  1981. 

[7]  Y.Y.  Lau  k  L.R.  Barnett,  “Theory  of  a  low  harmonic  field  gyrotron  (Gyromag- 
netron),”  Int.  J.  Infrared  Millimeter  Waves,  vol.  3,  pp.  619-744,  1982. 

[8]  P.S.  Rha,  L.R.  Barnett,  J.M.  Baird,  k  R.W.  Grow,  in  Int.  Electron  Device  Meet., 
pp.  525-538,  1985. 

[9]  R.W.  Grow  k  U.A.  Shrivastava,  “Impedance  calculations  for  traveling  wave  gy- 
rotrons  operating  at  harmonics  of  the  cyclotron  frequency  in  magnetron- type  cir¬ 
cuits  operating  at  the  PI  mode,”  Int.  Electron  Devices  Meet.,  pp.  384-387,  1982. 

[10]  H.S.  Uhm,  C.M.  Kim,  k  W.  Namkung,  “Theory  of  cusptron  microwave  tubes,” 
Phys.  Fluids,  vol.  27,  pp.  488-498,  1984. 

[llj  K.R.  Chu  k  D.  Dialetis,  “Theory  of  harmonic  gyrotron  oscillator  with  slotted 
resonant  structure,”  Int.  J.  Infrared  Millimeter  Waves,  vol.  5,  pp.  37-56,  1984. 

[12]  K.R.  Chu  k  D.  Dialetis,  “Kinetic  theory  of  harmonic  gyrotron  oscillator  with 
slotted  resonant  structure,”  Infrared  Millimeter  Waves,  vol.  13,  pp.  45-74,  1985. 

[13]  J.Y.  Choe  k  W.  Namkung,  “Experimental  results  of  cusptron  microwave  tube 
study,”  IEEE  Trans.  Nuclear  Sci.,  vol.  NS-32,  pp.  2882-2884,  1985. 

[14]  U.A.  Shrivastava,  R.W.  Grow,  P.S.  Rha,  J.M.  Baird,  k  L.R.  Barnett,  “Threshold 
power  transfer  for  the  gyrotron  and  peniotron  oscillators  operating  at  the  harmonic 
cyclotron  frequencies  using  coaxial  electron  beam-circuit  configurations,”  Int.  J. 
Electron.,  vol.  61,  pp.  33-r59,  1986. 

[15]  P.  Vitello,  Eleventh  International  Conf.  Infrared  Millimeter  Waves,  Tirrenia,  Pisa, 
Italy,  pp.  46-48  (1986). 

[16]  A.M.  Ferendeci  k  C.C.  Han,  “Theory  of  high-harmonic  rectangular  gyrotron  for 
TEmn  modes,”  IEEE  Trans.  Electron  Devices,  vol.  ED-31,  pp.  1212-1218,  1984. 

[17]  C.C.  Han  k  A.M.  Ferendeci,  “Nonlinear  analysis  of  a  high-harmonic  rectangular 
gyrotron,”  Ini.  J.  Electron.,  vol.  57,  pp.  1055-1063,  1984. 

[18]  A.M.  Ferendeci  k  C.C.  Han,  “Linear  analysis  on  an  axially  grooved  rectangular 
gyrotron  for  harmonic  operation,”  Int.  J.  Infrared  Millimeter  Waves,  vol.  6,  pp. 
1267-1282,  1985. 


74 


k 


[19]  Q.F.  Li  k.  J.L.  Hirshfield,  "Gyrotron  traveling  wave  aimplifier  with  out-ridged 
waveguide,”  Int.  J.  Infrared  Millimeter  Waves,  vol.  7,  pp.  71-98, 1986. 

[20]  J.  D.  Jackson,  Classical  Electrodynamics,  (John  Wiley  k  Sons,  Inc.,  New  York, 
1962). 

[21]  G.  B.  Collins,  Microwave  Magnetrons,  (McGraw-Hill,  New  York,  1948). 

[22]  M.  Abramowitz  k  A.  Stegun,  Handbook  of  Mathematical  Functions,  (Dover  Pub¬ 
lications,  Inc.,  New  York,  1972). 

[23]  G.F.  Brand,  “A  gyrotron  frequency  detuning  equation,”  Int.  J.  Infrared  Millimeter 
Waves,  vol.  4,  pp.  919-931,  1983. 

[24]  P.  Vitello  k  K.  Ko,  “Mode  competition  in  the  gyro-peniotron  oscillator,”  IEEE 
Trans.  Plasma  Sei.,  vol.  PS-13,  pp.  454-463,  1985. 

[25]  P.  Vitello,  W.H.  Miner,  k  A.T.  Drobot,  “Theory  and  numerical  modeling  of  a  com¬ 
pact,  low-field,  high-frequency  gyrotron,”  IEEE  Trans.  Microwave  Theory  Tech., 
vol.  MTT-32,  pp.  373-386,  1984. 


i 


75 


Figure  Captions 

Figure  la.  Cross-section  of  the  slotted  cylindrical  gyrotron. 

Figure  lb.  Cross-section  of  the  slotted  rectangular  gyrotron. 

Figure  2.  Guiding  center  coordinate  system. 

Figure  3.  Es  RF  field  for  a  ridged  cylindrical  cavity  with  Rg  =  1.4,  =  ir/2N,  iV  =  6. 

Figure  4.  Variation  of  the  start  oscillation  condition  for  the  third  harmonic  electron  cy¬ 
clotron  maser  interaction  (dashed  curves)  and  second  haurmonic  peniotron  inter¬ 
action  (solid  curves)  as  a  function  of  Ro  for  the  TEsn  mode.  The  beam  axial 
velocity  is  =0.1. 

Figure  5.  RF  field  for  a  ridged  rectangular  cavity  with  L*  =  2.186,  Ro  =  1.198,  w  = 
0.05466,  and  iVT  =  3. 

Figure  6.  TEoai  mode  start  oscillation  condition  for  the  rectangular  cavity  for  the  electron 
beam  with  72  keV,  and  =  0.28.  The  solid  curve  gives  QPt  for  the  ridged  cavity, 
while  the  dashed  curve  gives  QPb  for  the  smooth  walled  cavity. 
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